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ABSTRACT 

The s-process should occur in all but the lower mass progenitor stars of plan- 
etary nebulae, and this should be reflected in the chemical composition of the gas 
which is expelled to create the current planetary nebula shell. Weak forbidden 
emission lines are expected from several s-process elements in these shells, and 
have been searched for and in some cases detected in previous investigations. 
Here we extend these studies by combining very high signal-to-noise echelle spec- 
tra of a sample of PNe with a critical analysis of the identification of the emission 
lines of Z>30 ions. Emission lines of Br, Kr, Xe, Rb, Ba, and Pb are detected 
with a reasonable degree of certainty in at least some of the objects studied here, 
and we also tentatively identify lines from Te and I, each in one object. The 
strengths of these lines indicate enhancement of s-process elements in the central 
star progenitors, and we determine the abundances of Br, Kr, and Xe, elements 
for which atomic data relevant for abundance determination have recently be- 
come available. As representative elements of the "light" and "heavy" s-process 
peaks Kr and Xe exhibit similar enhancements over solar values, suggesting that 
PNe progenitors experience substantial neutron exposure. 

Subject headings: ISM: abundances — nuclear reactions, nucleosynthesis, abun- 
dances — planetary nebulae: general 



1. Introduction 

As remnants of stars that have evolved through the asymptotic giant branch (AGB) 
phase, most planetary nebulae (PNe) are believed to consist of material that has undergone 
nuclear processing in the precursor star via the s-process. The analysis of nebular emission 
from elements that have experienced nucleosynthesis in the parent star provides valuable 
information for stellar models. However, the detection of emission lines from ions enhanced 
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by the s-process has been hampered by their weakness and by uncertainties in the atomic 
data needed for the analysis, including line wavelengths. 



An initial attack on this problem was made a decade ago by iPequignot &: Baluteau 



( 11994J ) (hereafter PB94), who obtained a deep optical spectrum of the high-ionization PN 
NGC 7027. Using the best atomic data available for the energy levels of the more prominent 
ionization stages of elements in the fourth and fifth rows of the periodic table, they identified 
a number of post-Fe peak emission lines in the nebula. From the observed line intensities 
they concluded that the elements Kr and Xe, the latter normally a predominantly r-process 
element in stars of solar metalicity, were enhanced in NGC 7027 by factors of ~20 relative 
to their initial formation abundances, presumably roughly solar. They also detected Ba II 
and [Br III] emission at intensities indicating that Ba could be enhanced whereas Br might 
be depleted relative to solar values, subject to uncertainty due to poorly known excitation 
cross sections. 

Dinerstein and collaborators have subsequently pursued the study of s-process abun- 
dances in PNe through su r veys to detect IR fine-stru cture nebular emission from post-Fe 
peak ions (jDinerstein!l200ll ; ISterling fc Dinerstein! 120041 ) and from far-UV resonance-line ab- 
sorption by s- process elements in t he intervening nebular shell seen in FUSE spectra of their 
central stars (ISterling el al.ll2002t ISterling &: Dinerstein! 120031 ) . They detected IR emission 
from Se and Kr in roughly 50% of their sample PNe from which abundances of up to 10 
times solar were derived for these elements. The FUSE spectra revealed Ge III absorption in 
five PNe for which Ge enhancements of 3-10 times solar were deduced, clear evidence that 
central star AGB progenitors are major sites of s-process element production. 

As part of our on-going program to detect and identify the weakest lines in PNe at 
visible wavelengths down to levels significantly below that of the continuum, viz., < 10~ 5 the 
intensity of H/3, we have obtained very high signal-to-noise spectra at high spectral resolution 
of some of the higher surface brightness nebulae. Some fraction of the weaker lines observed 
are likely to originate from post-F e peak elements and therefore we have used the automatic 
line identification routine EMILI (jSharpee et al.ll2003l ). supplemented by recent energy level 
data for heavier elements, to assist in the search for such lines. Since the publication of PB94 , 
new calculations have be en made for spontaneous emission coeffici ents (IBiemont et al.lll995l ) 
and collision strengths (jSchoningj 119971 ; ISchoning fc Butlerl Il998l ) for fourth and fifth row 
elemental ion transitions. In this study, we use these newer atomic parameters to compute 
ionic and overall elemental abundances for Br, Kr, and Xe in order to make comparisons 
with their abundances in H II regions and the sun, and to derive s-process enrichment factors 
relevant to the study of s-process nucleosynthesis in the progenitor stars. 



2. Observations and Data Reduction 



Our present sample of objects consists of four PNe (IC 2501, IC 4191, NGC 2440, and 
NGC 7027) which satisfy the most important criteria for the detection of weak emission 
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lines in having (a) high surface brightnesses, (b) low expansion velocities, and (c) except for 
NGC 7027, relatively small internal dust extinction. The latter two criteria produce sharper 
lines and higher peak intensities relative to the continuum. 

We obtained spectra of the PNe IC 2501, IC 4191, and NGC 2440 during two observing 
runs of two nights each in early 2003 using the Las Campanas Observatory (LCO) Baade 
6.5m telescope with the MIKE echelle spectrograph. Similar instrumentation setups were 
used during the two runs. MIKE is a dual-beam spectrograph which simultaneously measures 
separate red and blue spectra, giving useful data over the continuous wavelength ranges 3280- 
4700 A and 4590-7580 A at resolutions (A/AA) of 28,000 and 22,000, respectively, for the 1 
arcsec slit width we used. 

A series of long and short exposures were taken of each nebula on one or more nights, 
typically adding up to a period of 2-3 hours at a time spent observing each object. Because 
MIKE is used without an image rotator at a Nasmyth focus, the orientation of the slit with 
respect to each nebula rotated by a large amount during these series of exposures. The central 
star was placed along a line perpendicular to the spectrograph slit so that the slit center 
was roughly midway between the central star and outer edges of the nebula, and then the 
telescope was guided to keep the star in that same position as seen on the acquisition/guide 
TV. The result was that during the course of the observations the slit swept out an arc about 
the central star, so that the final spectrum integrates over an area of the nebula that is much 
larger than the 1x5 arcsec slit. The total integration times in the combined long exposures 
were 630 min for IC 2501, 450 min for IC 4191, and 330 min for NGC 2440. A journal of 
the observations is given in Table [IJ and includes the approximate position of the center of 
the slit with respect to each of the central stars at the start of the series of exposures. 

The optics of the MIKE spectrograph introduce strong distortion into the image formed 
on the detector, so the projected emission lines have significantly different tilts as a function 
of their position upon it. When extracting spectra of objects such as these PNe which fill 
the slit, the tilts introduce unacceptable smearing in the wavelength direction (typically 3 
pixels over the 40 pixel slit length) unless they are corrected for during the extraction of the 
one-dimensional spectra from the two-dimensional echelle image. We wrote our own set of 
auxiliary FORTRAN programs to calibrate the tilts and perform the proper extraction. This 
procedure was tested on the emission lines in the comparison lamp spectrum. Compared to 
the emission line from a single pixel at the slit center, lines summed over the full slit length 
came out at the same pixel location in the dispersion direction and were broadened by 2 
percent on average. We are thus confident that the effects of these tilts are negligible in 
the extracted PNe spectra. These tilt-corrected spectra were then fed into the same suite of 
IRAF-based reduction programs used with the NGC 7027 data as described next. 

NGC 7027 was observed on 4 nights in June 2002 with the Mayall 4m Telescope at Kitt 
Peak National Observatory using the Cassegrain echelle spectrograph. Because we were 
interested in detecting faint 40-50 km s _1 FWHM emission lines over the widest possible 
wavelength range, rather than in detailed measurements of the line profiles, we used the 
short UV camera with the 79.1 grooves mm -1 echelle grating, cross-dispersing grating 226-1 
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in first order, and a GG-475 order separating filter. This gave full wavelength coverage from 
4600 to 9200 A at ~20 km s~ x FWHM (A/AA = 15,000) resolution with our 2 arcsec slit 
width, with partial coverage out to ~9900 A. On each night we offset to the same position 
with the slit at PA 145° and centered 0.5 arcsec N and 3.5 arcsec W of the central star, which 
is the brightest part of the PN shell in Ha emission. Most of the observing time was spent 
obtaining sequences of 1200s exposures which added up to a total of 680 min of integration 
time. We also took a number of shorter exposures 30, 60, 120, and 300s in length to measure 
bright emission lines that were saturated on the longer exposures. 

The data for all four PNe were reduced using standard procedures with IRAF-based 
reduction packages in the same manner that has been described in detail in our discussion of 
comparable sp ectra of the PN IC 4 18 that we obtained previously with the CTIO 4m echelle 



spectrograph f Sharpee et al. 2004h. A corre ction was made for the presence of Rowland 



ghosts near strong lines ([Baldwin et al.l 120001 ). To correct for flexure and temperature drift 
effects in the spectrographs, comparison lamp spectra were taken at roughly one hour inter- 
vals, and the wavelength calibration for each PN spectrum was made using the comparison 
lamp taken nearest to it in time. The wavelength fits indicate a 1-cr uncertainty in the wave- 
length scale that varies over each echelle order, and is of the order 4-6 km s _1 for all of our 
spectra. The spectr a were flux ca l ibrate d using observations of several spectrophotometric 
standard stars from lHamuy et al.l ( 19941 ) . 



The final extracted flux and wavelength calibrated spectra of all the PNe were used 
for our line identification and analysis, and are available in FITS format as an electronic 
supplement to the present article. 



3. Emission Line Selection 

The line selection procedure is initiated by fitting the continuum of each echelle order 
of the final spectrum with a smooth function that is pegged to the observed continuum 
at approximately ten wavelength intervals distributed over each order. We developed an 
automated procedure that selects what should be reliable continuum points for the fitting 
by selecting wavelength regions with a paucity of lines. Because bad pixels, artifacts caused 
by scattered light in the spectrograph due to strong emission lines, and line blends can 
cause poor fits, the fits were reviewed manually for each order and adjusted as necessary 
to insure that the continuum representation was valid. An example of the output of the 
automatic continuum fitting algorithm is shown in Figure [H where a typical fit prior to 
manual adjustment is displayed. Once the proper continuum level is established, each order 
is sampled pixel by pixel to detect emission features, which are defined as regions where the 
observed flux exceeds the continuum flux by la or more over a wavelength interval equal to 
or greater than that of the resolution of the spectrograph. 

All putative emission features were examined individually by eye and compared with 
their appearance on the original two-dimensional echelle images to establish their reality. We 
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have found that all features with fluxes greater than Ylo of that of the continuum, i.e., with 
S/N> 12, are clearly visible on the echelle images. Since scattered light features usually 
trail across multiple orders one can generally distinguish such artifacts from real lines by 
visually examining the images. The most uncertain aspect of the line detection process is 
distinguishing multiple line blends in real emission features. The expansion velocities of PNe 
are such that intrinsic line widths do vary by factors of 2-3 with the level of ionization, so 
it can be difficult to discriminate between one broad line vs. two or more closely spaced 
narrow lines. Figure [T] gives an example of the features within a selected wavelength region 
in one of our PNe spectra that have been designated by our software as emission lines. 

For detection of the weakest nebular lines, where the distinction between a continuum 
noise spike and a real feature can be difficult to establish with certainty, we have inter- 
compared the spectra of those PNe that have similar levels of ionization but different radial 
velocities. Noise spikes do not generally produce features having the width of the instru- 
mental resolution, and instrument artifacts tend to occur on the same place on the detector, 
which is at different rest wavelengths in the PNe spectra because of their differing radial 
velocities. Since spectra were obtained on different instruments, (MIKE and the 4-meter 
KPNO echelle), all significant scattered light ghosts could be detected through comparisons 
between the PNe spectra obtained with MIKE and the NGC 7027 spectrum. Our final line 
lists do contain a few weak features having S/N<7 which are present in at least two of the 
objects. 

The telluric nightglow emission spectrum was also sampled by these spectra, as sky 
subtraction in these extended objects was deemed impractical, particularly in regards to the 
subtraction adding significant additional noise. However, most nightglow lines were distin- 
guishable on the two-dimensional spectra by their uniform intensity and characteristic shape 
and size, namely those of the imaging slit. The positions and intensities of prospective night- 



glow l ines were compa r ed to those listed in the telluric feature atlases of lOsterbrock et al. 



( 119961 ) and iHanuschikl (120031 ). and likely matches removed from the nebular list except in 
those cases where their blending with stronger nebular lines was deemed only a minor con- 
taminant. For lines considered as candidate s-process ion transitions, nightglow s ources were 



given additional scrutiny through a comparison with the original spectra o f the lHanuschik 



( 120031 ) atlas, the likely identifications of the atl as lines (ICosbv et al.l 120061) . time-averaged 



observed intensities of their constituent systems ( Cosby fc Slanger 20061 ). and in some cases 



model spectra of those systems constructed with the molecular simulation code DIATOM 



(http://www-mpl.sri.com/software/DIATOM/DIATOM.html). 



Following the selection of sets of likely nebular lines in each of the PNe, their wavelengths 
and intensities were determined by single or multiple, in the case of blended features, Gaus- 
sian fitting of their profiles and immediate underlying continua. Simple summing was used 
for the most irregular and ill-defined line profiles with intensity-weighted centroids utilized 
as wavelengths. Emission lines appearing in multiple orders were then collated and their 
attributes averaged together. Emission line intensities and wavelengths from short and long 
duration exposure spectra of each PN were then normalized to a particular fiducial frame. 
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Line wavelengths were shifted to the nebular rest frame through a comparison with 
either the Balmer (MIKE PNe) or Paschen (NGC 7027) sequence l aboratory wave lengths. 
Line intensities were de-reddened with the Galactic extinction law of iHowarthl (119831 ). utiliz- 
ing an iterative process involving the magnitude of either the Balmer or Paschen jump, and 
the Balmer or He II 5-n sequence decrements to establish electron temperatures and c(H/3) 
logarithmic extincti on at B.0 values. The c(H/3) values for each PN are in good agreement 
with those listed in ICahn. Kaler. fc Stanghellinil (jl992j). Errors in the final line intensities 
in all PNe spectra, as determined from the formal errors to the profile and continua fits, are 
similar to those found for NGC 7027: 41% for I<10 5 I(H/3), 20% for I=10- 5 -10~ 4 I(H/3), 
11% for I=10~ 4 -10 -3 , and >6% for I>10 -3 I(H/3) on average. This is independent of any 
errors arising from the reddening correction. For weak lines, where the greatest contributor 
to uncertainty is the indeterminate level of the true continuum, the formal errors in the 
fit probably understate the actual uncertainty in the intensity. From random inspection of 
several lines at the 10~ 5 H/3 intensity level, it is likely that some lines may have uncertain- 
ties ranging upwards to 100% in regions where the continuum level is rapidly undulating, 
complicated by scattered light artifacts from adjacent orders, or affected by strong telluric 
absorption. Final geocentric offsets and c(H/3) values are presented in Tabled 

Figure [2] presents a histogram showing the fraction of lines in each of the nebulae ob- 
served with MIKE for different intensity levels relative to B./3 that we have determined to be 
real and which also appear in at least one of the other two PNe (the NGC 7027 spectrum is 
omitted here due to its different wavelength coverage and spectral resolution). The strongest 
emission lines are all detected in each of the nebulae. Even at intensities down to 10~ 5 H/3, 
where S/N~10 typically, roughly 50% of the weak features identified in the individual PNe 
also appear in one of the other nebulae, suggesting that they are true nebular lines. The 
present spectra are among the deepest emission spectra taken of nebulae, revealing some of 
the weakest lines yet observed. This can be seen in Figure [31 where the cumulative number 
of lines exceeding a given flux level relative to H/? is shown for several of the most extensive 
PNe spectral studies published in the recent literature. 



4. Plasma Diagnostics and Abundances 



Electron te mperatures and densiti es are presented in Table [31 The IRAF nebular pack- 
age task temden ( iShaw fc Dufourlll995l ) was used to equate relative intensities of collisionally 
excited diagnostic lines to densities and temperatures by matching a diagnostic for each from 
ions of similar ionization energy, and solving for self-consistent values. Each ion was modeled 
by a five or greater level atom for this purpose, with spontaneous emission coeffic ients and 
collisi on strengths for the five lowest levels drawn primarily from the compilation of lMendoza 
( 119831 ). Although these are not the default values currently distributed with nebular, this 
atomic data set yielded good agreement in temperature and density among diagnostics in 
our previous analysis of IC 418, and were utilized in the o r iginal nebular release. Spon- 
taneous emission coefficients from iFroese Fischer &: Tachievl (120041 ) and collision strengths 
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from IWilson fc Belli (120021 ) were used for N I and CI II respectively. For the remaining en- 
ergy levels and for the ions CI IV and K V, the atomic data used in the most recent release 
of nebular was utilized. The departure of the [K V] density diagnostic from other density 
diagnostic values, particularly in NGC 2440 and IC 4191, may indicate errors in this atomic 
data. However, s-process elemental abundance derived later using this diagnostic (Kr +4 /H + 
and Xe +5 /H + ) were completely insensitive to the utilized density. Errors in diagnostic values 
were determined by selecting the extrema values from computations at all combinations of 
diagnostic line ratios plus and minus their (1-cr) uncertainties, including an error estimate 
for the reddening correction. Indeterminate error limits, such as occurred when a ratio value 
exceeded the asymptotic limit or where the paired diagnostics failed to converge at a partic- 
ular ratio value, are listed without a value. Balmer (MIKE PN e) and Paschen (NG C 7027) 
jump temperatures were calculated in the manner described by IZhang et al.l (120041 1 . and as 
reported by them for many PNe, they are lower than those temperatures determined from 
the collisionally excited lines. 

Ionic abundances derived from both collisionally excited and nominal radiative recom- 
bination lines are presented in Table HI The IRAF nebular package tasks ionic and abun- 
dance were used to make computations from the collisionally excited lines listed in Table [3] 
using temperatures and densities derived from the diagnostics with the closest ionization 
potential to each ion. Where no diagnostics where clearly appropriate, averaged values for 
temperature and density where used. Uncertainties were computed in the same manner 
as the diagnostic values, by calculation of the abundance at every combination of tem- 
perature, density, and line intensity ratio plus or minus their (1-cr) uncertainties (where 
available) and selecting the extrema values. Uncertainties were not calculated for abun- 
dance computed with averaged diagnostics values. For the recombination lines, effective 
recombination coefficients were combined with line intensities to make abundance determi- 



nations, in the manner described by ISharpee et al.l (120041 ). with temperature and densities 
again drawn fro m diagnostics ne arest in ionization potential to each ion. For He + /H + , 
the coefficients of ISmita (119961 ) or lBenjamin. Skillman. &: Smita ( 119991 ) were used to deter- 
mine an average abundance from the A4923, A5876, and A6678 lines, with corrections for 
collisi onal excitation (Case A triplets, Case B for single ts) taken from iKingdon fc Ferland 
( 119951 ). while coefficients from [Storey fc Hummerl (119951 ) were used to calculate an average 
He +2 /H + abundances from various He II lines. To determine C +2 , N +2 , Q +2 , and Ne +2 abun 



danc e s relative to H + , recent effective recombination coefficients (IStoreylll994l ; lKisielius et al. 



19981 ; iDavey. Storey, fc Kisieliusll2000l ; iKisielius fc Storey! 120021 ) for the strongest observable 
multiplets were used to calculate abundances. As seen in Table 0J +2 /H + and Ne +2 /H + 
abundances deduced from the collisionally excited lines are systematically lower than those 
derived from radiative recombin ation lines, as is generally the case in PNe (jLiu 
Robertson- Tessi fc Garnettll2005l ). 



2004 
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Emission Line Identification 



The emission line identification code EMILI (ISharpee et al.l 120031 ) was used to make 
the majority of emission line identifications. EMILI creates models of the ionization-energy 
dependent velocity field and ionization level of a PN or H II region from user-supplied 
empirical data. These models are used by EMILI to select from a large atomic transition 
database all transitions within five times an observed line's wavelength measurement error, 
and to compute relative intensities for emission lines corresponding to those transitions. 
Those transitions predicted to produce emission line intensities within 3 orders of magnitude 
of the highest value among all transitions initially selected are then subjected to a test for 
the presence of lines from the same LS-coupled multiplets (but not all lines from the same 
upper level) at expected wavelengths and relative intensities. Potential identifications are 
then assigned a numeric likelihood parameter based on their wavelength agreement with the 
observed line, strength of the predicted emission line, and results of the multiplet check, and 
are ranked and presented. 

Many emission lines from Z>30 elemental ions were observed in NGC 7027 by PB94. 
To place these lines on an equal footing for identification purposes with those arising from 
more abundant lighter elements, EMILI used the Atomic Line List v2.05 of P. van Hoofl as 
its reference database, which extends to Z=36 (Kr). The latest experimental determinations 
available in the literature for ground electron configuration energy levels of Z>36 ions were 
then added to this database. Transition wavelengths for the mostly optically forbidden 
transitions among these levels were constructed by differencing the level energies. Sources 
for all Z>30 ion energy levels and associated atomic data used in subsequent analysis are 
provided in Table All ions with certain or probable line identifications in NGC 7027, 
as rated by PB94, most of those with more tenuous identifications, and a handful of other 
ions they suggested to be worthy of future consideration at a higher spectral resolution, were 
incorporated into the database. However, ions with level uncertainties perceived or explicitly 
stated to be greater than 1.0 cm -1 in their source literature were excluded. 

EMILI was run against each nebula's set of observed emission line wavelengths and 
intensities. Electron temperature and density values were provided to EMILI by the results 
of standard plasma diagnostics, derived from strong collisionally-excited lines with certain 
identifications. Construction of the empirical kinetic and ionization models were also drawn 
from those same lines. The EMILI reference elemental abundances were set to the solar 



values of iLodderd (120031 ). EMILI was also run against the IC 418 line list of ISharpee et al. 
(120031 ). to determine if some of its remaining unidentified lines could be identified with an 
expanded set of transitions, and to also act as a foil for the mostly higher excitation PNe 
considered in the present sample. The line identification lists produced by EMILI were then 
visually inspected, and the EMILI-preferred assignments were compared to identifications in 
the literature for those same lines in PNe and H II region spectra of similar depth, spectral 



1 http://www. pa. uky.edu/peter/newpage/ 
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resolution, and level of ionization. The entries in our final line list are in many cases taken 
from among the highest-ranked EMILI suggested identifications. 

A segment of a final line list is given in Table El the full line list is available in the 
electronic version of this manuscript. This table lists all observed lines with their character- 
istics, viz., observed wavelength in the nebula rest frame, identifications, full widths at half 
maxima, S/N, observed and reddening corrected intensities relative to H/5, and their most 
likely identifications. Lines with uncertain identifications are denoted by a likely blends 
with a "bl" or "ns" if blended with a telluric emission feature, and uncertain features are 
marked with a "?" . Only the perceived strongest component of a blend is listed in the table. 
Lines without obvious identifications are listed here without an identification. 



6. Z> 30 Line Identifications 

PB94 detected 25 emission lines from several Z>30 ions in their optical spectra of 
NGC 7027, including Se, Br, Kr, Rb, Sr, and Y from the fourth row of the periodic table, 
Xe and Ba from the fifth row, and Pb from the sixth. Eighteen of these detections were 
considered "certain" or "probable", and 7 considered "possible". Lines from 13 additional 
Z>30 ions were also either tentatively identified or proposed as future targets for spectra 
with greater spectral resolution and better S/N. 

Given the depth and high-resolution of the spectra considered here, confirmation of 
the PB94 identifications in multiple PNe was sought, as were additional lines belonging 
to other Z>30 ions. The use of EMILI allows prospective Z>30 transitions and weaker 
transitions of more abundant lighter elements to be treated equally for purposes of emission 
line identification. Collision strength and spontaneous emission coefficient calculations for 
these transitions allow accurate predictions to be made of the relative intensities of lines 
arising from the same ion, allowing identifications of their lines to be made with more 
confidence. 



The line lists for the present PNe sample and IC 418 (ISharpee et al.ll2003l ) were searched 
for the strongest expected Z>30 ion transitions within their observed bandpasses. These 
primarily were the 3 Pi,2 _1 D2 nebular transitions of ions with 4p 2 ' 4 and 5p 2,4 valance electrons, 
the 4 Sg , 2 - 2 Dg / 2 5 , 2 nebular transitions of ions with 4p 3 and 5p 3 valence electrons, and the 



po 2po 

1/2,3/2 r1 



"l/2 3/2 fine-structure transitions for 4p, 4p , 5p, 5p , and 6p valance electron ions. 
Fine structure 3 P- 3 P transitions, when accessible in the visible, were also included. For the 
cases of Kr IV, Xe IV, and Br IV, auroral 2 ^>y 2 5/2~ 2 ^i/2 3/2 anc ^ trans-auroral ^s/2" 2 ^ > i/2 3/2 
transitions were also considered. The permitted resonance lines of Ba II and Sr II were also 
searched for. Initially, all wavelength coincidences of 1A or less between an observed line and 
a transition wavelength were considered possible Z>30 lines, regardless of the identifications 
recommended by EMILI for those lines. 



Tables \7\- fTUl present excerpts from the EMILI output for observed lines believed to rep- 
resent the best cases for a Z>30 transition as the actual identification for an observed line in 
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at least one of the PNe in which that line appeared. For each observed line in each PN, the 
identifications of highest rank, as specified by the "Identification Index" or IDI, the EMILI 
figure-of-merit for quality of an identification (5 or less is considered a quality identification), 
and their multiplet search statistics (numbers expected/observed, hereafter multiplet statis- 
tics) are included. The IDI value is followed by a letter A— >D if among the top four highest 
ranked transitions, with "A" being the highest. Additional identifications drawn from the 
literature and from the terrestrial nightglow (prefaced with "SKY" ) are also included. Iden- 
tifications that did not yield predicted line intensities within three orders of magnitude of 
the strongest value among all identifications, or were outside the 5-cr search radius, do not 
have a calculated IDI value. The IDI value/rank is sometimes followed by a symbol, an as- 
terisk indicating the most likely single identification, "bl" indicating components of a likely 
blend, and ":" indicating indeterminate alternate identifications. The reddening corrected 
intensity of the feature attributable to the putative s-process identification, if corrected for 
any blending as described in succeeding sections (due to other identifications denoted with 
a "bl"), is presented in italics below its originally observed value. The limiting intensity, 
determined from the local minimum flux of line detection considered certain (S/N=7), or 
from imposition of artifical lines of S/N=7 at that wavelength, is presented in brackets for 
the case of s-process transitions without corresponding observed features. An "OUT" label 
indicates features residing outside the observed bandpass of a spectrum. Finally, features of 
dubious reality are denoted by a "?" following the observed wavelength. 

Figure @] depicts continuum-subtracted spectra, with wavelengths shifted to the nebular 
rest frame, in the vicinity of those Kr, Xe, and Br lines for which a positive identification 
was made in at least one PN, and which were potentially observable in all four PNe of the 
present sample and in IC 418. Comments on identifications pertaining to individual Z>30 
elemental ion lines follow. 



6.1. Kr Line Identifications 

PB94 noted that [Kr III] 4p 4 3 P 2 -4p 4 X D 2 A6826.70 has long been observed m various 
novae and PNe, but has seldom been identified as such. They identified [Kr III] as a strong 
line in NGC 7027, blended with weak C I (V21) A6828.12 and He I 3s 3 S-16p 3 P° A6827.88 
lines. 

As seen in Table [7] and Figure [5], the higher resolution of our present PNe spectra cleanly 
separates both the He I and C I line from the putative [Kr III] line, except for NGC 7027, 
where He I contributes minimally to its red shoulder. However, Figure |5] shows that the 
R branch head of the telluric nightglow OH Meinel 7-3 band is a serious contaminant in 
this region. In NGC 7027, a broad O VI A10 32 Raman scattering line at A6829.16, seen 



previously in NGC 7027 by IZhang et al.l (120051 ). is also observed. The OH contribution was 



represented by a model of the band normalized in intensity to the uncontaminated Ri(1.5) 
line at A6834.01. The He I 3s 3 S-15p 3 P° A6855.91 line profile was shifted and scaled by a 
factor of 0.83 to represent He I A6827.88 (Smits et al. 1991; Case B, for the most appropriate 
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grid point: T e =10 4 K, n e =10 4 cm -3 ). The profile of the companion O VI Raman scattering 
line observed at A7088, scaled upwards by a factor of four to yield the best fit with the 
red wing of the putative [Kr IV] feature, represented the A6829.16 line. These line profiles 
and the telluric model were subtracted from the continuum normalized spectra to yield the 
residual features shown in Figure [51 

For all PNe except NGC 2440, where the observed line appears to be entirely comprised 
of the OH line, there appears to be a substantial residual feature that is coincident, but 
except for IC 4191, slightly to the red of the predicted wavelength for [Kr III] A6826.70. 
The EMILI results for the corresponding observed lines suggest that besides the C I and 
He I transitions mentioned earlier, the [Fe IV] A6826.50 feature is probably the only other 
sensible alternative identification. However, given the absence of other multiplet members 
that should be observ ed with this lin e and the relative observed strengths of the likely 
strongest [Fe IV] lines (jRodrfguezll2003l ) . this is not a likely identification except perhaps for 
IC 4191. Unfortu nately, the compani o n [Kr I II] 3 Pi- x D 2 line at A9902.3, predicted to be ~18 
times weaker by IBiemont fe Hansen! (jl986bl ). is only potentially observable in NGC 7027, 
where a line listed at A9903.55 is much too strong and most likely associated with C II 4f 
2 F°-5g 2 G A9903.67. Nevertheless, we conclude that the [Kr III] A6826.70 identification is 
relatively certain in IC 418, IC 2501 and NGC 7027, uncertain in IC 4191, and that line is 
not observed in NGC 2440. 

Examination of Table [7] shows that both transitions of the [Kr IV] 4 S°- 2 D° nebular mul- 
tiplet AA5346.02,5867.74 are clearly identified in every spectrum, appearing as the primary 
EMILI ranked identifications. For A5346.02, the [Fe II] A5347.65 identification is clearly 
too far away and the S II (V38) A5345.71 dielectronic transition is unlikely given the ab- 
sence of other multip let members. The Al II A 5867.8 identification for A5867.74, given by 
Baldwin et al.l (120001 ) and ISharpee et al.l (12003), is clearly sup erseded. With a theoretical 
intensity ratio I(A5346.02)/I(5867.74) of 0.65 (jSchoninjdll997l ). assuming electron densities 



well below the critical density values of ~18-1.3xl0 7 cm -3 for the 2 D^ 2 and 2 D^ 2 levels 
respectively, as justified from the values in Table the observed A5346.02 intensity appears 
to be somewhat too high for IC 4191 and IC 418. This suggests the possible blending of 
[Kr III] A5346.02 with another transition, perhaps C III (V13.01) A5345.85, although poor 
multiplet statistics and the weakness of the strongest C III lines in IC 418 cast doubt on this 
possibility for that PN. However, the other three PNe have ratios of 0.76 (NGC 2440), 0.72 
(IC 2501), and 0.73 (NGC 7027), which are slightly higher but consistent with the expected 
value within the combined measurement errors of the line intensities. 

The auroral [Kr IV] 2 D° /25/2 - 2 P° /2 AA6107.8,6798.4 identifications are less straightfor- 
ward, except for NGC 7027 where both are considered certain according to both PB94 and 
the present EMILI results. In NGC 2440, the observed l ine corresponding to the A 6798.4 
transition, which should be 2.5 times weaker than A6107.8 (IBiemont fe Hansenlll986bl ). is ac- 
tually stronger, while in IC 2501 the stronger A6107.8 line is not present at all. As suggested 
by both PB94 and the EMILI results, C II (V14) A6798.10 may be responsible for either 
the total intensity in IC 2501 or excess intensity in NGC 2440. While PB94 note that C II 
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(V14) is of dielectronic recombination origin and its multiplet components should not have 
relative intensities expected due to LS-coupling rules, an examination of the intensities of 
all lines from this multiplet appearing in each PN suggests that the intensities do appear to 
roughly follow those rules with the exception of what would be the strongest line at A6783.91. 
Therefore, as an approximation, a LS-coupled relative strength of A6798.10 to A6791.47 of 
0.16, instead of the one-to-one ratio with A6812.28 used by PB94, was utilized to correct the 
putative [Kr IV] A6798.4 feature in both NGC 2440 and NGC 7027, as A6812.28 appears it- 
self appears in only one of the PNe spectra (it would be the weakest feature in the multiplet 
if LS-coupling rules held). Subsequent abundance analysis also suggests that the [Kr IV] 
A6107.8 line itself may be too strong relative to its nebular counterparts in NGC 2440 and 
IC 4191, admitting [Fe II] A6107.28 as an alternate identification that appears more likely 
in the latter PN. In summary, both auroral transition can only be identified comfortably in 
NGC 7027. 

The third magnetic dipole transition of the [Kr IV] auroral multiplet 2 Dy 2 - 2 P°/ 2 A7131.3 
is not definitively detected in any spectrum including NGC 7027. However, in the NGC 7027 
spectrum, a feature is observed at A7131.65 with intensity 7.7 x 10~ 5 I(H/3), comparable to 
A6107.8. It originally was assumed to be an under-corrected Rowland ghost arising from 
the nearby saturated [Ar III] 3 P 2 - 1 D 2 A7135.80 line, but might instead be associated with 
A7131.3. The electric quadrupole line [Kr IV] 2 D5/ 2 ~ 2 Pi/ 2 A8091.0 is forty times weaker and 
therefore undetectable. 

PB94 claim a detection of [Kr V] 'Pi^Da A6256.06 as a blend with C II (V10.03) 
A6257.18 at their instrumental resolution. In our spectra, as is shown in Figure [5j C II 
A6257.18 is resolvable from the putative [Kr V] feature in all spectra except NGC 7027, 
where it s lightly contr i butes to its red wing. However, in IC 2501 and NGC 7027, as also 



noted by IZhang et all (120051 ) . another C II line, dielectronic C II (V38.03) A6256.52, is 
believed to be a significant contaminant in some PNe based upon the strength of the nearby 
A6250.76 line from the same multiplet and upon favorable EMILI assessments. The profiles 
of both C II lines were approximated and subtracted from all the PNe spectra in Figure O 
where the C II (V10.03) A6259.56 line profile scaled downward by a factor of 0.56 represents 
A6257.18, and the C II A6250.76 profile, also scaled downward by a factor of 0.56 according 
to LS-coupling statistics, represents A6256.52. The contribution of the interloping telluric 
OH 9-3 band was accounted for through subtraction of a model of the band normalized in 
intensity to the nearby Qi(2.5) A6265.21 line. 

The resulting residual plots show a distinct line just to the red of the predicted wave- 
length of [Kr V] A6256.06 in NGC 2440, IC 2501, and NGC 7027, that we believe can be 
definitely identified as such in all three PNe. For IC 4191, the residual flux, 3.0 x 10~ 6 I(H/3), 
is probably too low to be an actual line, and the entirety of the original profile is ascribed 
here to C II (V38.03) A6256.52, particularly since both the original line and A6250.76 share 
the same broad profile in the spectrum. In IC 418, the putative [Kr V] profile is completely 
removed by the combination of the nightglow and C II (V38.03) A6256.52 model profiles, as 
is appropriate given the nebula's low excitation. 
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The other nebular line, [Kr V] 3 P2- 1 T>2 A8243.39, accessible only in the bandpasses of 
the NGC 7027 and IC 418 spectra, sits amidst the head of the Paschen series, rendering it 
difficult to disentangle at the resolution of our NGC 7027 spectrum. In the IC 418 spec- 
trum, a comparable observed line is indistinguishable from other nearby lines in the Paschen 
series and is clearly identifiable as H I 3-43 A8243.39, with the Kr V identification unlikely 
due to the PN's low ionization level. Returning to NGC 7027, PB94 identified this tran- 
sition at A8242.7 as a blend with N I (V2) A8242.39 and O III 5g G 2 [9/2]°-6h H 2 [11/2] 
A8244.10, compromised by telluric absorption. In our NGC 7027 spectrum, [Kr V] A8243.39 
is tentatively identified as the blue peak at A8244.33 of a resolved two-line blend with H I 
3-42 A8245.64, and which appears significantly affected by telluric absorption. The alter- 
nate identification for the line, N I A8242.39, is clearly resolved here as a separate feature. 
Because this observed line has a measured intensity ratio with A6256.06 roughly comparable 
to what would be expected if both lines were due to [Kr V], I(A6256.06)/I(A8243.39)=1.1 
( jBiemont fc Hansen! Il986al ). it is believed that O III A8244.10 contributes only a minor 
amount to the line, although the presence of the telluric absorption complicates matters. To 
summarize, both [Kr V] lines appear to be present with about the expected intensity ratio 
in the spectrum of NGC 7027, but both are absent in the spectrum of IC 418. For the cases 
of IC 2501 and NGC 2440, [Kr V] A6256.06 is probably present. 

The spectra were also examined for evidence of other auroral and trans-auroral tran- 
sitions of Kr ions. Only the auroral line [Kr V] 1 D 2 - 1 So A5131.78 satisfied the initial ±lA 
screening criterion and also appeared in the EMILI output as a possible identification for 
an observed line at ~5131.0A in all PNe. The case for this identification is lessened by 
the expected weakness of this line compared to other [Kr V] lines that were not definitively 
identified, and particularly because of its high observed intensity of 1.3 x 10~ 3 I(H/5) relative 
to the ~ 1 x 10~ 4 I(H/3) for AA6256.06,8243.39 in NGC 7027, the nebula that might be 
expected to exhibit the best evidence for this line given the strength of the other confirmed 
Kr lines. Instead either O I 3p 3 P-8d 3 D° A5 13 1.25, which was the highest ranked line in all 
but one nebula (IC 418) and that has a strength comparable to other members of the same 
sequence, or C III 5g 3 G-7h 3 H° A5130.83 appear the more likely identifications. 



6.2. Xe Line Identifications 

A number of Xe ion transition identifications are also considered probable in our spectra. 
The EMILI statistics for Xe identifications are affected by th e low solar Xe abundance, 



which is an order of magnitude lower than Kr (ILoddersI 120031 ) . This contributes to low 
predicted relative emission line intensities and consequently higher or absent IDI values for 
its identifications as compared to those computed for weak transitions from more abundant 
elements. As such, every appearance of a Xe ion transition as a candidate identification for an 
observed line in an EMILI list, regardless of its IDI value, was given serious consideration as 
it signaled that alternative identifications from more abundant elements did not predominate 
despite the advantage arising from greater abundances. Table [8] lists the EMILI statistics 
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for lines judged most likely to correspond to Xe ion transitions, while Figure H] depicts the 
regions around the most likely observed transitions. 

The identification of [Xe III] A5846.77 in NGC 7027, was given by PB94 to an 

excess intensity in the He II 5-31 A5846.66 line. This excess was searched for in the present 
PNe sample through the subtraction of the He II A5837.06 profiles, shifted to the A5846.66 
line position and scaled assuming I(A5837.06)/I(A5846.66)=0.92 (Storey and Hummer 1995; 
Case B, for the most appropriate grid point: T e =10 4 K, n e =10 4 cm -3 ) as is shown in FigureEJ 
The intensity of He II A5837.06 was not corrected for the presence of C III 7h 3 H°-18i 3 I 
A5836.70 because the line used for this correction by PB94, C III 7i 3 I-18k 3 K° A5841.2, was 
either not present or was an identification of low rank (IDI=6, ranked 5th) for a corresponding 
line in NGC 7027. As seen in Figure El distinct residual profiles are present after subtraction 
of He II A5846.66 for both NGC 7027 and IC 418 (the latter having negligible He II), 
and at the correct wavelength for [Xe III] A5846.77. A bizarrely shaped profile is seen for 
IC 2501, and no profiles are seen for either NGC 2440 or IC 4191 after subtraction. While 
the corresponding observed lines are probably He II A5846.66 in the latter two PNe, IC 2501 
has very weak He II lines, suggesting that the residual profile is attributable to something 
else, tentatively [Xe III] A5846.77. The [Xe III] A5846.77 identification, while not appearing 
in the EMILI lists for NGC 7027 or IC 2501, is considered a better choice than the other 
top ranked identification, [Fe II] A5847.32, which has poorer wavelength agreement and 
multiplet statistics in all cases. While the [Xe III] 3 Pi- 1 D2 transition at 1.37/i is unavailable 
for confirmation, we believe that the [Xe III] A5846.77 is definitely present in NGC 7027 and 
IC 418, but is only tentatively identified in IC 2501. 

The reality of the observed features potentially corresponding to the [Xe IV] 4 S^ 2 3 , 2 - 
2 Dgy 2 AA5709.2,7535.4 transitions is marginal except in NGC 7027. However, their identifica- 
tions as the Xe lines, should they be actual emission lines, is more certain. For A5709.2, the 
widely observed N II (V3) A5710.77 line is detected as a separate line in all of the spectra. 
Alternative identifications such as [Fe I] a5D3-a5Pi A5708.97 either have too many missing 
multiplet members, except for IC 2501 where the ratio of potentially observed to total num- 
ber of multiplet lines expected is marginally better, or as in the case of the Fe II] A5709.04 
intercombination line, unlikely given that permitted Fe II lines were not anywhere defini- 
tively identified. For A7535.2, the reality of the observed lines is less uncertain, including 
IC 418 where a previously unreported line has been uncovered by more tho rough exami 



nation of its spectra. The Fe II] (V87) A7534.82 identification proposed by iHyung et al. 



fcOOlt ) is unlikely for the same reasons as Fe II] A5 709.04 for [Xe IV] A5 709.2 The N II 5f G 
2 [7/2l 4 -10d 4 F° A7535.10 identification, selected bv lEsteban et all tooi ) and lPeimbert et al. 



( 120041 ) , also appears unlikely given the relative strengths of other known permitted N II lines 
in the spectra. It should be noted that the EMILI statistics for the A7535.10 identifications 
in the MIKE PNe sample would probably have been better if the line did not appear near 
the end of the last spectral order where the wavelength calibration is the poorest. Assuming 
an electric density well below t he critical regime (~ 10 7 cm -3 ), the expected intensity ratio 
is I(A5709.2)/I(A7535.4)=0.64 JSchoning fe Butlerlll998h . The observed range of values that 
we observe, 0.33-1.19, suggests measurement errors on par with the average values for such 
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weak lines. The strongest examples in NGC 7027 do yield the best agreeing ratio (0.82). As 
such, we conclude that both [Xe IV] lines may be present in all PNe spectra with varying 
degrees of certainty. 

One of the stated goals for future spectroscopy of NGC 7027 and other PNe given by 
PB94 was the detection of [Xe V] transitions, none of which they detected with any certainty. 
In the present spectra, no observed line passed the ±lA initial selection criterion for the 
3 P 12 - 1 D 2 AA5228.8,6998.7 lines, nor did these transitions appear as possible EMILI IDs for 
any observed line. This is despite the improved sensitivity and resolution which should have 
enhanced the chances of detecting A5228.8, and which should have easily separated A6998.7 
from O I (V21) A7002.10, the contaminant noted by PB94. Nevertheless, the detection 
of these lines in our spectra would remain problematic since for A5228.8 extensive flaring 
from [O III] A5006.84 in an adjacent order leads to numerous ghosts in its vicinity, while 
for A6998.7 telluric absorption in the tail end of the Fraunhofer A band complicates its 
observability. 

EMILI did suggest that the fine-structure transition [Xe V] 3 P - 3 P2 A7076.8 as a possible 
identification for observed lines in two PNe, IC 4191 and NGC 7027, PNe with appropriate 
excitation levels for the appearance of a line of this ion. This is the primary EMILI identi- 
fication in IC 4191. A third occurrence of the corresponding observed line in NGC 2440 is 
clearly associated with a nightglow line from the O2 b 1 S+-X 3 E~ 3-2 band upon comparison 
with a simulation of that band. The competing identification C I (V26.01) A7076.48 has 
poor multiplet statistics, while [Fe III] A7078.10 and [Ni II] A7078.04 are feasible, but have 
poor wavelength agreement. The Ca I identifications are also unlikely given that they arise 
from energy levels of large energy, where lines that would follow from cascades from these 
levels are not observed. Two remaining obstacles to an identification with [Xe V] A7076.8 are 
the lack of the nebular 3 P- 1 D lines and the low branching ratio for this electric quadrupole 
transition (~0.008 with respect to 3 Pi- 3 P 2 at 2.07/x). Nevertheless, despite this transition's 
expected weakness, the lack of viable alternate identifications suggests the corresponding 
observed line may be at least tentatively identified in NGC 7027, and probably identified in 
IC 4191, and both occurrences yield reasonable abundance values in subsequent analysis. 

The fine structure transition [Xe VI] 2 Pi/2~ 2 Ps/2 A6408.9 was identified with certainty 
in NGC 7027 by PB94. In the present NGC 7027 spectrum, A6408.9 is probably associated 
with a weak but clearly separable observed line on the red wing of He II 5-15 A6406.38. The 
alternate primary EMILI identification of [Fe III] A6408.50 was not considered likely, as an 
inspection of the NGC 7027 spectrum for other lines from all low energy [Fe III] multiplets 
did not show a significant number of matches to warrant strong consideration as a likely 
identification. [Fe III] A6408.50 also has a comparatively high excitation energy (6.25 eV). 
A C IV 9-17 A6408.70 identification is similarly downgraded by its high excitation. The 
reality of corresponding observed lines in NGC 2440 and IC 4191 is questionable, although 
for lack of suitable alternate identifications, uncertain identifications of [Xe VI] A6408.90 are 
retained for both. 

As with Kr, an inspection for auroral and trans-auroral transitions of various Xe ions was 
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also undertaken. There were coincidences within ±lA between the [Xe III] 3 Pi- 1 So A3799.96, 
[Xe III] ^a-^o A5260.53, [Xe IV] 4 S° /2 - 2 P° /2 A3565.8, [Xe IV] 2 D° /2 - 2 P° /2 A4466.5, [Xe IV] 
2 D° /2 - 2 P° /2 A5511.5, [Xe IV] 2 D° /2 - 2 P° /2 A6768.9, and [Xe V] W-^o A6225.3 lines. How- 
ever, except for two instances, the identifications did not appear in the EMILI list for those 
particular lines, and most had more reasonable and higher ranked, identifications: C II (V30) 
A5259.66,.76 for example instead of the auroral [Xe IV] A5260.53 line. In many cases the line 
only appeared in one of the PNe spectra analyzed here, and could be rejected due to the 
strength or absence of the nebular transitions from the same ions. 



6.3. Br Line Identifications 



In PB94 two transitions of Br were identified: [Br III] 4 S3^ 2 - 2 Dgy 2 A6131.0 with certainty, 
and [Br IV] 2 D° /2 - 2 P° /2 A7385.1 as possible. A second nebular transition, [Br III] 4 Sjj /a - 2 D|? /2 
A6556.4, was lost in a blend amidst the [N II] A6548.04+Ha+[N II] A6583.46 complex. 



The source for the Br III levels in the Atomic Line List v2.05 is iMoord (119581 ). while 
the source used by PB94 is stated to be the experimental levels listed in Table IV of 
Biemont &: Hansen! (Il986al ) from unpublished work of Y.N. Joshi and Th.A.M. van Kleef, 
and provided by private communication with van Kleef. However, a comparison of the Ritz- 
determined wavelengths from both sets of levels with those listed by PB94 shows that the 
Joshi and van Kleef 2 D° term energy levels were used for the nebular transition wavelengths, 
while the IMoord (119581 ) levels were used for auroral transit i on wa velengths. The substantial 
difference in the 2 D3^ 2 level energy, 15042 cm -1 for IMoord (119581 ) and 15248 cm -1 for Joshi 



and van Kleef, leads to a difference in the 4 S3 / / 2 - 2 D^ 2 transition air wavelength of 6646.3 A 
versus 6556.4 A, res pective l y, with a lesser difference for the other transition, 6132.9 A versus 
6131. OA. Since the IMoord ( 19581 ) levels used by the Atomic Line List v2.05 and therefore 
by EMILI h ave stated uncertainties of 0.63 cm -1 , while those for Joshi and van Kleef in 
Table IV of iBiemont fc Hansen! (11986a! ) are listed to a precision 1 cm -1 , nearly the same 
level of uncertainty, it is difficult to know which levels are more accurate. Therefore, EMILI 
was run using both sets of energy levels, and the 4 S°- 2 D° lines were searched for at both sets 
of resultant wavelengths. 



For the transitions wavelengths generated from the IMoord (119581 ) energy levels, A6132.9 
and 6646.3, no observed line was detected in any spectrum meeting the initial ±lA selection 
criterion. However, for the 4 S3/ 2 - 2 D^ 2 transition wavelength generated from those energy 
levels used by PB94, A6131.0, there do appear to be observed lines in 4 of the 5 PNe spectra 
near the A6130.4 line that PB94 identified with [Br III]. As seen in Table El while [Br III] 
A6131.0 is the highest ranked EMILI identification in only one PN. 

However, other identifications for the putative [Br III] A6131.0 line are not compelling. 
PB94 cite their observed line as being a blend with C III 7h 1,3 H°-16g 1,3 G A6130.30, which 
is the highest EMILI-ranked transition in each spectrum in which the putative Br III line is 
observed, except for IC 2501. However, the companion line C III 7h 1,3 H°-16i 1,3 I A6126.30, 
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claimed by PB94 to have an equal intensity to A6130.30, is not present in either the IC 4191 
and IC 418 spectra. Well-known lower excitation C III lines are of negligible intensity in 
IC 418, while in IC 2501 C III A6126.30 has poor wavelength agreement with its corresponding 
observed line, and is not the primary EMILI identification for that line. Only in NGC 7027 
is C III A6126.30 a primary ID, suggesting that only in its spectrum is the equally intense 
C III A6130.30 likely to be present and accountable for at least a portion of the putative 
[Br III] A6131.0 line. As such for NGC 7027, the A6126.30 line intensity is subtracted from 
the observed line at 6130. 32A, while in the remaining spectra the corresponding observed 
lines are ascribed wholly to [Br III] A6131.0. Possible contamination due to flaring from 
an adjacent order in the NGC 7027 spectrum was not accounted for. The [Ni VI] A6130.40 
identification, the second ranked identification in many of the spectra, is of high excitation 
(8.3 eV upper level), too high an ionization for IC 418, and appears only to be ranked highly 
due to a favorable coincidence in wavelength. 

Spectroscopic confirmation for this identification is sought in the possible presence of 
the [Br III] 4 S3/ 2 - 2 Dg\ 2 A6556.4 transition, in the PNe spectra at expected intensities similar 
to that of A6131.0. As mentioned, detection of this line is difficult given its proximity to 
the saturated Ha line and its associated ghosts in its immediate proximity. Nevertheless, 
inspection of the line lists and spectra (Figure HJ) does show distinct co-aligned features in the 
IC 2501, IC 4191, and NGC 7027 spectra near this wavelength. While initially assumed to 
be ghosts, the fact that these features can be found in spectra from two different instruments 
and appear invariant to the Ha intensity, suggest they may be unrelated to Ha. That real 
lines can be observed between [N II] A6548 and Ha is demonstrated by the presence of the 
strong telluric OH Pi (3.5) 6-1 A6553.62 line just to the blue of the suspected [Br III] feature. 
EMILI lists the [Br III] A6556.4 identification among possible choices in all PNe in which 
the line appears, with the A6131.0 identification satisfying the multiplet search in IC 2501 
and IC 4191. Alternate EMILI-favored identifications for these features correspond to O II, 
and N II permitted and core-excited transitions between levels of primary quantum number 
4-6, and are doubtful as they are of comparatively high excitation with respect to better- 
known 3-3 or 3-4 transitions in these PNe that have equal or lesser intensities. They appear 
here only due to their wavelength agreement with the observed lines. The highest ranked 
Fe II A6555.94 is discounted once again by the lack of other permitted Fe transitions in any 
spectra. 

The observed intensity ratios of the two putative [Br III] lines, I(A6131.0)/I(A6556.4), 
ranging from 0.18-0.45, do conflict with theoretical expectation as follows. While specific 
collision strengths are unavailable at present for the Br III states, the ratio of the strengths 
relevant to these transi tions should still be roughly proportional to their respec tive upper 



level statistical weights (jPequignot fe Baluteaulll994l ; lOsterbrock fe Ferlandll2006l ). weighted 



by a Boltzman factor respecting the difference in level energies, and non-negligible fine- 
structure emission between the 2 D° levels, at sub-critical electron densities. Employing the 
IRAF nebular package task ionic to solve for the relative populations of the 2 D° levels at 
temperature derived from the diagnostic thought to be the most appropriate, from [Ar III], 
and using appropriate Kr IV collision strengths as a proxy, scaled as discussed in Section 



-18- 



the I(A6131.0)/I(A6556.4) ratio is expected to be 0.72-0.83. This is close to the 0.67 value 
expected from the ratio of the collision strengths alone, and two to three times larger than 
the ratio observed in the spectra. At these levels the [Br III] A6131 line should have been 
observable in NGC 2440, as the predicted intensity exceeds the likely detection limit for 
features in its vicinity. Therefore, we consider [Br III] to have been only tentatively detected 
in NGC 2440. However, at least one [Br III] line does appear to be present in the remaining 
four PNe. They all have lines detected at the same wavelengths, although shifted 20-30 km 
s _1 to the blue of the wavelengths we have adopted for the [Br III] lines, at approximately the 
same relative intensities, and lack satisfying alternate identifications. The ratio of observed 
intensities (or in the case of IC 418 the upper limit), while not exact, are within a factor of 
2-3 of those expected. 



The identifications of the nebular [Br III] lines suggests that the level energies of iMoore 

(119581 ) for at least the 2 D° term levels may be in error. PB94 also identified an observed 



line at 7384. 3A as possibly auroral 2 Dy 2 - 2 P^ 2 A7385.2, but no such line appears in any 
of the PNe spectra examined here. Ho wever, s i nce P B94 appears to have used 2 p°- 2 D° 



transition wavelengths derived from the iMoord (119581 ) energy levels, this might not be a 



surprise, although no lines in any PN were discovered at the correspon ding wavelength for 
the 2 D? /2~ 2 P3/2 transition (7483.1 A) computed from the levels listed by lBiemont fc Hansen 



(11986ah either. 



The nebular [Br IV] 'P^D 2 A7368.0 line was detected by PB94 at 7366.0A, affected 
by telluric absorption, and blended with C IV 10-21 A7363.9 and O III 4s 3 P?-4p 3 Pi 
A7365.35. In our NGC 7027 spectrum, a line appearing at 7367.62A as a small protrusion 
above the local continuum level interpolated between telluric absorption features in both 
spectral orders covering that wavelength, is tentatively identified as [Br IV] A7368.0. The 
alternate candidates listed by PB94, dielectronic C II 3p' 2 D 5/ 2-3d' 2 P° /2 A7377.00 is too 
far away, while a O II A7367.68 identification is doubtful given its high upper level energy. 
The C V 7p 3 P°-8d 3 D A7367.60 transition is a possible alternate identification, but only 
one other C V transition in the spectrum, C V 6gh 1,3 G,H°-7hi 1,3 H°,I A4944.50 with an 
IDI value of 3, is a top ranked EMILI IDI; others are of lesser rank or do not appear in the 
EMILI lists for the corresponding observed lines. 

The other nebular transition, [Br IV] 3 P2- 1 D2 A9450.5 was not detected by PB94, even 



though the line should be of the same intensity (IBiemont fc Hansen! Il986af ). Comparison 
between the IC 418 and NGC 7027 spectra shows that something is filling in the telluric 
absorption feature in the latter at 9450. 85A, with an estimated intensity close to [Br IV] 
A7368.10, although a ghost feature at this wavelength arising from scattered light within 
the spectrograph cannot be ruled out based on the proximity of other similar features. If 
the feature is real, however, the alternate identification of Fe I A9450.95 is not compelling. 
Therefore, both transitions are identified in NGC 7027, although tentatively since they ap- 
pear in only one PN spectra, may be attributable to misinterpretation of the local continuum 
level complicated by telluric absorption of varying degree, and the latter observed line may 
be spurious. No auroral or trans-auroral features of [Br IV] were identified in any spectrum. 
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6.4. Other Z>30 Line Identifications 

The spectra were searched for lines originating from other Z>30 ions. Among the nu- 
merous coincidences between observed lines and transition wavelengths at the ±1A level, 
those listed in Table [TU1 were judged to be the most likely to correspond to real lines describ- 
able by a Z>30 identification in at least one of the PNe. Given the low abundances of these 
ions, only transitions among the lowest-lying levels were expected to be observable. 

PB94 declared as certain the identification of a line at 5758. 7A in NGC 7027 as [Rb IV] 
3 P 2 - 1 D 2 A5759.44. Corresponding lines are detected in 4 of 5 PNe here. However, at even 
the highest instrumental resolution, this transition would blend with He II 5-47 A5759.74. In 
NGC 2440 and IC 4191 there is contamination from a flare or ghost at the position of He II 5- 
47 A5759.74 in one spectral order. Using only the uncontaminated order, and inspecting the 
nearby He II 5-n sequence suggests that in IC 4191 the He II A5759.74 line shows significant 
excess, while in NGC 2440 there is better agreement with recombination theory and the line 
is likely He II. The NGC 7027 spectra, taken with a different instrument, does not show 
any contamination, but does show a similar excess in He II 5-47 A5759.74. Thus, the line 
intensities in IC 4191 and NGC 7 027 were corrected by 1( 5-47) /I(5-46)=0. 95 derived from 
the ratio of the their emissivities ( Storey fe Hummer 19951 ). IC 418 has no detectable He II 



lines, enhancing the [Rb IV] A5759.44 identification, but the weakness and irregularity of 
its profile cast doubt upon its reality as a line, and it is only tentatively identified here. 
The alternate EMILI recommended identification, [Fe II] A5759.30, was not considered likely 
given its high excitation energy. Since the expecte d intensity ratio of 3 P 2 - 1 D2 A5759.44 
to sPi^Da 9008.75 is 17 faiemont fe Hansen! [l986bh . the apparition of A9008.75 in IC 418 



and NGC 7027 is probably likely due to He I 3d 3 D-10p 3 P° A9009.23,.26, or is unreal. In 
summary, the [Rb IV] A5759.55 line is identified in IC 4191, NGC 7027, and tentatively in 
IC 418. 

Another Rb line detection considered probable by PB94 in NGC 7027 is [Rb V] 4 S° /2 - 
2 Dgy 2 A5363.6, which they identified at 5364. 2A. Although it does not appear in the EMILI 
list for the corresponding lines in our spectra, we believe that its identification in at least 
NGC 7027 is viable given the poor multiplet statistics of [Ni IV] A5363.35, where none of 
the four other multiplet lines are clearly present. The O II 4f F 2 [4]° / , 2 -4d / 2 F 7 / 2 A5363.80 
identification is an interesting alternative. Well-known O II dielectronic doublet lines of 
multiplets V15, V16, and V36, such as 3s' 2 D 5/2 -3p' 2 F° /2 (V15) A4590.97 are all clearly 
present in all of our PNe spectra, except NGC 7027 where they either outside the bandpass 
or not optimally placed for detection, at intensities ~ 10~ 4 I(H/3 ) . The se lines are present 



at similar intensities in the NGC 7027 spectrum of IZhang et all (120051 ). Lines from 3d-4f 
transitions, particularly 3d 2 D 5/2 -4f F 2 [4]° /2 (V92a) A4609.4 arising from the lower level 
of the O II A5363.80 transition are also favorably identified at roughly the same intensity. 
Therefore it is not out of the realm of possibility that this line is evidence of the partial 
feeding of 2 [4]£/ 2 level. However, while EMILI did not perform a multiplet check on this line 
(due to a non LS-coupled lower level), two other transitions at 5361. 74A and 5375. 57A from 
the same upper term and ending on this level are not present in any of the spectra. Without 
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spontaneous emission coefficients it is difficult to judge a potential branching ratio for these 
transitions. Therefore, it is believed that [Rb IV] A5363.60 might be tentatively identified 
in IC 2501 and IC 4191 as well. Unfortunately the matching 4 S° /2 - 2 D° /2 A4742.40 line is 
disguised by a flare or ghost in the MIKE spectra, and is not optimally placed in NGC 7027 
to allow a spectroscopic confirmation, nor observed in IC 2501 and IC 4191. 

Numerous lines of various Sr ions were searched for, including the strong (I/I(H/3)=3 x 
10~ 4 ) [Sr IV] 2 P3/ 2 ~ 2 Pi/2 A10276.9 fine-structure line named as a possible identification by 
PB94 in NGC 7027, [Sr VI] 4 S° /2 - 2 D° /2 A4249.2 and 2 Dg /0 - 2 P° / A5434.3 identifie d as pos- 



sible and tentative respectively by PB94 and also identified by IZhang et al.l (120041 ). and the 



Sr II resonance lines at AA4 077.71,4215.52 and the [Sr II] 2 Si /2 - 2 D 5/ 2,3 /2 AA6738.39,6868.17 



lines observed in rj Carinae (IZethson et al.ll200ll ). While, there were some instances of wave- 



length coincidences with observed lines, the magnitude of the wavelength differences and 
existence of plausible alternate identifications, such as [Fe II] 4249.08 for [Sr VI] A4249.2, 
did not warrant a claim of identification. The same was true for Y and Zr, where the [Y V] 
2 P° /2 - 2 P° /2 A8023.6 and [Zr VII] 3 P 2 - 3 P A7961.4 fine-structure lines , listed as possible and 
tentative identifications by PB94 respectively, did not have any wavelength coincidence with 
an observed line in any PNe within the initial ±lA screening criteria. 

Turning to the fifth row of the periodic table, we believe that an observed and previously 
unidentified line in IC 418 may correspond to [Te III] 3 Pi- 1 D 2 A7933.3, which would be the 
first optical identification of a Te III line in a PN. The He I 3p 3 P°-ns 3 S sequence of lines, of 
which the alternate identifications He I AA7932.36,.41 are members, does not become clearly 
evident in this spectrum until 3p-10s at 8632. 76A and 8632. 83A. The wavelength coverage of 
our spectra does not extend out to the possibly stronger 3 P 2 - 1 D 2 A10876.0 transition, so this 
line cannot be used to check the Te III identification. Similarly, the detection of a possible 
[I III] 4 S° /2 - 2 D° /2 A8536.6 line in NGC 7027, with a less than compelling Cr II permitted 
line as an alternate primary EMILI identification, cannot be confirmed through observation 
of its companion A Sy 2 - 2 D^ 2 A6708.7 transition. The putative A6708.7 line appears to be 
better explained by a combination of [Mn II] A6709.93 and possibly [Cr V] A6709.8, although 
the latter line has poor multiplet statistics (2/0). Nevertheless, the [I III] A8536.6 transition 
has a definite IDI that is second ranked for its corresponding observed line, so we count this 
as a tentative detection for NGC 7027. 

PB94 claims the detection of four transitions belonging to permitted multiplets of Ba II, 
with three detections considered certain, and one c onsidered possible . Some of these same 



transitions h a ve als o been identified in emission by (IZhang et al.ll2005l ) in NGC 7027 and by 



Hobbs et al.l (120041 ) in the compact H II region within the Red Rectangle. In the present 
sample it is believed that one of the transitions considered certain by PB94, Ba II 5d 2 Ds/ 2 - 
6p 2 P^/ 2 A6141.71, may correspond to an observed line in three of five PNe. Given the low 
solar abundance and high condensation temperature (1455 K; Lodders et al. 2003) of Ba, the 
detection of these transitions might initially be considered unlikely. However, as originally 
proposed by PB94, if sufficient gas-phase Ba + is available for collisional excitation, emission, 
and moderate self absorption, it is estimated that optical depths of 0.5 and 2.7 in the 6s 
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2 Si/ 2 -6p 2 ?3/ 9 A4554.03 resonance transition are sufficient to account for both its apparent 
non-detection, and for the observed intensity of the putative A6141.71 lines in NGC 2440 
and IC 4191 respectively. PB94 estimated an optical depth of 3 for the A4554.03 transition 
in NGC 7027. Under these scenarios the A6141.71 line is either among the strongest or is the 
strongest observable Ba II line, with its relative intensity with respect to other Ba II lines 
increasing with a larger optical depth in the A4554.03 line. 

The EMILI results do suggest some potentially viable alternative identifications other 
than Ba II for individual PN, but none that can satisfactorily account for the observed 
line in all three. O I 18d 3 D 2 -4p' 3 D3 A6141.75 arises from an auto-ionizing level, is too 
strong with respect to other permitted O I transitions further down the cascade chain, and 
is not accompanied by any other multiplet members in any PNe. Under the abundance 
and ionization model created by EMILI for NGC 2440 and IC 4191, Ne III 4p 5 P 3 -4d 5 D° 
A6141.48 did not produce an emission line with an intensity within three orders of magnitude 
of the strongest predicted intensity among all putative identifications. This was not the case 
in NGC 7027, where the line was among the strongest predicted lines, is expected to be the 
among the strongest in the multiplet, and where the identification is enhanced by favorable 
multiplet statistics. The [Ni III] 4s 3 F 2 -4s 3 Pi A6141.83, arises from a level of high excitation 
energy (9.9 eV), is probably not the strongest member of its multiplet as it does not originate 
from the level of the multiplet with the highest statistical weight, and other potentially 
stronger multiplet members are not observed in IC 2501 and NGC 7027. Yet in NGC 2440, 
the multiplet statistics are somewhat better, and the large intensity of the corresponding 
observed line and the Ba abundance derived from it is at odds with lower abundances derived 
for Kr and Xe in subsequent abundance analysis. The excellent agreement between the 
observed wavelengths in all three PNe and the Ba II A6141.71 wavelength, the expectation 
that A6141.71 is the strongest Ba II line, and the lack of a good alternate identification 
for IC 4191, suggests that Ba II A6141.71 warrants serious consideration as the correct 
identification in all cases. 

The fine-structure transition [Ba IV] 2 P3/ 2 ~ 2 Pi/ 2 A5696.6, another certain detection 
from PB94, is verified by its detection in our NGC 7027 spectrum after a correction for C III 
(V2) A5695.92 is made. The intensity attributable to A5696.6 was determined using the 
effective recombi nation coefficient for C I II A56 95.9 2 (a = 3.1 x 10~ 15 cm 3 s" 1 ) specified by 



PB94, those from iNussbaumer fc Storey! (119841 ) and lPequignot. Petitjean. &: Boissonl (1199ll ) 



for C III 5g 1,3 G-6h 1,3 H° A8196.61, and the A8196.61 line's observed intensity. Subtraction 
yielded a line amounting to 36% of the originally observed intensity. 

The only sixth-row elemental transition sought in these spectra was the fine-structure 
[Pb II] 2 Pi /2 - 2 P3 /2 A7099.80 line, which was identified with certainty by PB94. The transition 
appears in the EMILI list for IC 418 tied for the highest ranked transition for a previously 
unidentified line. For NGC 7027, the transition is not ranked, but the competing identifi- 
cations in its EMILI list are not convincing given their high excitation energies. Excellent 
wavelength agreement is seen in both cases, and both identifications are considered certain. 



In summary, of the 18 Z>30 elemental ion transitions considered certain or probable 
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by PB94 in NGC 7027, 15 are believed to be detected to various degrees of certainty in the 
present set of spectra. The final intensities of all Z>30 lines detected with any certainty 
within at least one PN are presented in Table HU with certain identification listed in bold, 
and tentative identifications in normal type. Conspicuous among those missing from PB94 
is [Se III] 3 Pi- 1 D 2 A8854.2. The nearby He I 3d 3 D-llp 3 P° A8854.14 transition does appear 
to show a large excess relative to that expected from other confirmed lines in the same series 
in both IC 418 and IC 70 27. However, the energy levels in the Atomic Line List v2.05, 
derived from iMoord (jl952j) and also utilized by PB94, have a listed uncertainty of 6.3 cm -1 
that exceeds the maximum 1 cm -1 tolerance allowed for inclusion of supplemental Z>36 
ions into EMILI. This leads to an error of up to 10A in the transition wavelength. Thus, 
the identification cannot be confirmed under the present degree of energy level uncertainty. 
Also missing are the Ba II transitions 6s 2 Si/ 2 -6p 2 P^ 2 A4554.03 and 5d 2 D 3 / 2 -6p 2 Pg , 2 1 , 2 
AA5853. 67,6496. 90. As discussed previously, all would be weaker in our spectra than the 
potentially observed A6141.71 under the assumption of a moderate optical depth in the 
resonance transition A4554.03 (and 6s 2 Si/ 2 -6p 2 P°/ 2 A4934.08 not observed by PB94). 



Zhang et al.l (120051 ) has claimed the identification of 5 additional transitions in their 



NGC 7027 spectra belonging to Z>30 elements that were not observed in our spectra. Two 
additional permitted Ba II lines, Ba II A4554.03 and Ba II A4934.08 are detected, both 
weaker than the A6141.71 line as would be expected for substantial optical depths of the 
resonance lines. They also identify auroral [Br III] A7385.2 and [Rb V] A5080.2, although the 
nebular [Br III] A6131.0 and [Rb V] AA4742.4, 536 3.6 lines, which mig ht be expected to be 
stronger than the auroral lines, are not identified. IZhang et al.l (120051 ) also identifies [Sr VI] 
A4249.2, which might be better attributable to [Fe II] A4249.08. 



6.5. Z>30 Line Identifications in H II Regions 



In comparing the spectra of the PNe with H II regions, we note that no post-Fe peak 
lines were identified by lGarcfa-Rojas et al.l (120041 ) in their deep VLT UVES echelle spectrum 
of the H II region NGC 3576. In Table [11] are included estimates of the upper limits to 
undetected Kr and Xe line intensities for this H II region from the faintest detectable lines 
that were identified at neighboring wavelengths or from artificial lines inserted at their wave- 
lengths that meet m inimal detection S/N statistics. In the spectrum of the Orion Nebula, 
Baldwin et al.l (120001 ) reported weak (I« 2 x 10 -5 H/?) unidentified features at 5867.8 A and 
6826.9 A that were originally unidentified, but can now be positively identified as [Kr IV] 
4 S° /2 - 2 D° /2 A5867.74 and [Kr III] 3 P 2 - 1 D 2 A6826.70 respectively. However, there are no other 
close coincidences between any lines listed in Tables l7HTm except for 5363. 34A which would 
correspond with [Rb V] 4 S^ 2 - 2 D^ 2 A5363.60, an unlikely identification given the unrealis- 
tically high degree of ionization (40 eV) for an H II region. The observed intensities and 
upper limits are also included in Table [TTJ 
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Z>30 Abundances 



The line intensities given in Table [IT] were used to compute abundances for Ba, Kr, 
Xe, and Br ions. The atomic data listed in Table [3] were formatted for inclusion in 5- 
level atom models for abundance analysis with the IRAF nebular task abundance, or for 
the cases of the fine structure lines [Xe VI] A6408.9 and [Ba IV] A5696.6, using a two-level 
atomic solution code. For the [Br III] and [Br IV] lines the relevant collision strengths 
have not yet been calculated. However, since these ions are isoelectronic with [Kr IV] and 
[Kr V], and because collision strengths fo r the same le vels along an isoelectronic sequence 
tend to vary with effective nuclear charge (jSeatonlll958l ). with some exceptions, the collision 
strengths of [Br III] and [Br IV] were assu med to be 25% smal l er tha n those for Kr. Because 
collision strengths were not calculated by ISchoning fc Butlerl (119981 ) for [Xe V] transitions, 
the [Kr V] collision strengths for the same transitions were utilized in the Xe +4 analysis. 
Collision strengths for excitation to the 3 P2 parent level of t he tentatively observed [Xe V] 
5p 2 3 P n - 3 P? A7076.8 line in other n p 2 ions such as Ne V (|Lennon fc Burke 1994 ). Ar V 
( IGalavis. Mendoza. fc Zeippenl Il995l ). and Kr V (ISchoningi 1 199 71 ). depart by 15% at most 
from their Kr V values at 10000 K. The temperatures and densities used for these analyses 
corresponded to those from diagnostics with the closest ionization potential. The [Ar III] 
(27.6 eV) temperature and [CI III] (23.8 eV) density diagnostic values were used for [Br III] 
(21.8 eV), [Kr III] (24.4 eV), [Xe III] (21.1 eV), and [Ba IV] (20.0 eV); [O III] (35.1 eV) 
temperature and [Ar IV] (40.1 eV) density for [Br IV] (36.0 eV), [Kr IV] (37.0 eV), and 
[Xe IV] (32.1 eV); [Ne III] (41.0 eV) temperature [Ar IV] (40.1 eV) density for [Xe IV] (46.0 
eV); [Ar V] (59.8 eV) temperature and [K V] (60.9 eV) density for [Kr V (52.5 eV)] and 
[Xe VI] (57.0 eV). The [O I] temperature and [N I] density were used for Ba II assuming 
collision excitation of Ba + as the source of Ba II A6141.71. Averaged diagnostic values were 
used when one of the diagnostics was unavailable for a particular ion. Uncertainties were 
computed in the same manner as the lighter elements, through permutation of intensity 
measurement and diagnostic value errors, where available, and selection of extrema values. 
For lines that were corrected for blending, 25% of the value of the correction was added to 
the measurement uncertainty. 

Table [T2l presents the results of the abundance determinations for individual lines of Ba, 
Br, Kr, and Xe ions. To compute overall elemental abundances for the latter three ions, argon 
has been selected as a benchmark element in addition to hydrogen. This was done because 
the ionization potentials of the noble gases Ar, Kr, and Xe, as well as Br, are very similar 
for their first three stages of ionization and therefore ionization corrections should not be 
large when making abundance comparisons among these elements. Additionally, the noble 
gases are almost completely non-reactive and have very low condensation temperatures, so 
corrections for gas phase abundances depleted by grain formation are insignificant for these 
elements. Only ionic abundances relative to hydrogen were computed for Ba. 



We convert from ionic to elemental abundances by making use of the similarity in 
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ionization potentials of the noble gases, so that 

Kr/Ar = (Kr +2 + Kr +3 )/(Ar +2 + Ar +3 ) and (Kr +2 + Kr +3 + Kr +4 )/(Ar+ 2 + Ar +3 + Ar +4 ) , 
Xe/Ar = (Xe +2 + Xe+ 3 ) /(Ar +2 + Ar+ 3 ) and (Xe+ 2 + Xe+ 3 + Xe +4 + Xe +5 )/(Ar+ 2 + Ar+ 3 + Ar+ 4 ) 

for H II regions and PNe judged to be of low excitation (IC 2501 and IC 418), and for high- 
excitation PNe (NGC 2440, IC 4191, NGC 7027), respectively. The Kr +4 /H+ abundance 
was included in the IC 2501 Kr/Ar ratio determination. For NGC 7027, the Br/Ar ratio 
was calculated the same way as the Xe/Ar ratio, with corrections for unobservable Br +4 
and Br +5 made assuming Br +4 ,Br +5 /Br=Xe +4 ,Xe +5 /Xe, appropriate because of their very 
similar ionization potentials. 

The resulting abundances ar e listed in the bottom rows of Table [12] relative to their 



solar values from iLodderd (120031 ). Three particular results are noteworthy. First, three 
out of five PNe show significant enhancements in Kr and Xe relative to solar values, while 
two others show similar, solar-like values. It is also interesting to note that IC 418 and 
NGC 7027, both considered young PNe, have the largest over-abundances, although they 
are of greatly different ionization classes. Meanwhile, the H II regions show only solar Kr 
abundances, indicative of unprocessed ISM gas. Secondly, the Kr and Xe abundances in 
all PNe all show enhancements of similar magnitude, as was seen in NGC 7027 by PB94. 
Finally, for NGC 7027, the Br abundance also shows a level of enhancement similar to that 
of Kr and Xe. Uncertainties remain regarding the adaptability of Kr collision strengths, 
modified as was done here, for the Br abundance calculations. However, to reduce the Br 
abundance to solar, the collisions strengths for Br +2 and Br +3 would have to both be four 
times greater than their counterparts in Kr, and this is much larger than the difference seen 
for analogous transitions for other p 2 and p 3 ion pairs (N and O, CI and Ar). In any event, 
the prevalence of probable Br line identifications suggests, independent of the actual atomic 
data, that significant Br does exists in most of the PN comprising our sample. 

It is of interest to compare the s-process abundances we derive for PNe with those 
obtained for evolved stars from analyses of their absorption spectra. With the exception of 
Rb, Ba, and Pb, the elements observed in emission in PNe are different from those normally 
observed in the spectra of late- type stars, so a direct comparison is not feasible. In fact, lines 
of the sam e s-process e lemen ts are not necessarily even observed in similar-type stars. For 



this reason iLuck et al.l (119911 ) defined two parameters, [Is] and [hs], that represent the mean 
abundance of elements associated with Sr and Ba, respectively, in what are called the "light" 
and "heavy" s-process peaks. They define the abundance indices [Is] and [hs] for an object as 
the mean logarithmic abundances relative to iron of (Y, Zr, & Sr) and (Ba, Nd, La, & Sm), 
respectively, compared to their solar mean abundances. All of these elements are produced 
by the s-process , although several of them have predominantly r-process contributions for 



solar- type stars (lArlandini et al.l 1 1 9991 ). 



Of the post-Fe peak elements that are observed in PNe, Kr belongs to the light s-process 
peak near Sr, and Xe is near Ba in the heavy s-process peak. Both are produced by the 
s-process, although Xe is predominantly an r-process element for stars of solar metalicity. 
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Using models to determine the appropriate correction factors, Kr and Xe can be incorporated 
into the [Is] and [hs] indices. However, in making comparisons of nebular abundances with 
those derived from stellar spectra Fe should not be used as the fiducial abundance for nebulae 
because its consistently strong depletion from th e gas phase due to grain formation c auses its 
true abundance in nebulae to be indeterminate (jShieldslll975l ; iPerinotto et al.lll999l ). Argon 
is a good surrogate for Fe because it does not suffer depletion in nebulae, and it is well 
observed in nebulae in multiple ionization stages and the relevant excitation cross sections 
are known. 

The extent to which s-process nucleosynthesis occurs in stars is determined by the total 
neutron exposure after the Third Dredge Up phase. Calculations show that neutron exposure 
is the primary factor that determines the ratio of the enhancements of the heavy to light 
s-process elements, i.e., [hs/ls], the "neutron-exposure-related-parameter". Low neutron 
exposures result in element production confined to a crowded region around the Sr peak, 
producing [hs/ls] < 0, whereas higher n e utron exposures significantly popula te the Ba peak, 
leading to [hs/ls] >0 (IBusso et al.lll995uBusso. Galling, fe Wasserburgj|l999l ). The fact that 
[Xe/Kr]wl for PNe, in spite of Xe not being produced by the s-process as much as Kr, 
is strongly suggestive that the progenitor AGB stars of PNe experience significant neutron 
exposure. Of equal significance, the fact that forbidden lines of post-Fe peak elements not 
observed in stars are detectable in nebulae causes H II regions and PNe to be of potentially 
great value in studying the nucleosynthesis of these elements. 



8. Summary 

In summary, very high signal-to-noise spectra of PNe do reveal lines from elements that 
are enhanced by the s-process, as was also found by other investigators. Because many 
post-Fe peak elements are refractory their gas phase abundances are not reliable indicators 
of the nucleosynthetic processes that have occurred in the progenitor stars. Fortunately, the 
noble gases Ar, Kr, and Xe are largely unaffected by molecular and grain formation, and the 
latter two elements are situated in the light and heavy s-process peaks, respectively. They 
also have easily excitable and observable lines in multiple ionization stages for which atomic 
cross sections are now available, and whose analyses should be straightforward. 

We find for a sample of five PNe that Kr and Xe abundances are enhanced over solar 
values by up to an order of magnitude, from both an analysis of their intensities relative to 
those of the fiducial element Ar, and also from a relative comparison of their line strengths in 
PNe vs. H II regions such as the Orion Nebula, where the lines are very weak in gas that is 
representative of the ISM composition. The similar enhancements of Kr and Xe in PNe are 
suggestive of large neutron exposures in the progenitor central stars. Further spectroscopy 
of PNe should reveal additional post-Fe peak emission lines whose analyses will contribute 
to a more complete picture of post-main sequence stellar evolution, while deep spectra of 
H II regions should lead to improved values of the ISM abundances of these elements 
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Table 1. Journal of Observations 



Object Telescope Dates(UT) Integration Times(sec) Slit Center 



NGC 7027 KPNO 4m 19-22 Jun 2002 34x1200s, 9x300s, 1x120s, 3.5" W, 0.5" N 

4x60s, 6x30s 

IC 2501 LCO 6.5m 12-13 Feb 2003 6x1800s, 1x300s, 4x60s 2" E 

IC 2501 LCO 6.5m 25-26 Feb 2003 15x1800s, 2x60s 2" E 

IC 4191 LCO 6.5m 12-13 Feb 2003 5x1800s, 1x1200s, Ix300s,lx60s 2" E 

IC4191 LCO 6.5m 25-26 Feb 2003 10xl800s,lx300s,3x60s 2" E 

NGC 2440 LCO 6.5m 12-13 Feb 2003 11 x 1800s, 1x300s, 1x60s 14" E 
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Table 2. Nebular Parameters 



Parameter 


NGC 2440 


IC 2501 


IC 4191 


NGC 7027 


V(geo) (km s _1 ) 


+66.7 


+23.2 


-40.1 


-7.1 


c(H/3) 


0.55 


0.55 


0.77 


1.17 
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Table 3. Electron Temperatures and Densities 



Diagnostic NGC 2440 IC 2501 IC 4191 NGC 7027 



Density (cm 3 ) 



[N I] A5198/A5200 


oouu_ 1300 


21000iii- 5 ooo 


12000rt b oo 


15000- 8000 


[S II] A6716/A6731 


3100™ 


11000 






[0 II] A3726/A3729 a 

b 


33001S8 


noooi°888 


9oool«888 




3700 


11000 






[CI III] A5517/A5537 


4700+ 1100 


85ooi?288 


12000™ 


47400+ 1900 


[Ar IV] A4711/A4741 c 


6200ti£8 


93ooi?288 


ii8ooi»88 




d 


5900liSgg 


7400™ 


11000™ 


49200™ 


e 


6300li«88 


99ooig88 


13000™ 




[K V] A4123/A4163 


43 00+ 19000 




41 000+ 17000 




Temperature (K) 










[0 I] (A6300+A6464)/A5577 


9400±300 


69001^8 


79001^88 


ii3ool 3 88 


[S II] (A6716+A6731)/(A4068+A4076) 


150001^88 


12000 






[0 II] (A3726+A3729)/(A7320+A7330) 


16000 


13000 






[N II] (A6548+A6583)/A5755 


126001™° 


10800l°°8 


12500±1100 




[Ar III] (A7136+A7751)/A5192 


mool™ 


94001^88 


130001*88 


12900±200 


[0 III] (A4959+A5007)/A4363 


147001^88 


95001^88 


9900±300 




[CI IV] (A7531+A8046)/A5323 


129001188° 


61001888 


86001^88 


13700±200 


[Ne III] (A3869+A3968)/A3343 


155001^88 


11000±200 


ii6ool 3 88 




[Ar V] (A6435+A7005)/A4625 


16400™ 




neoo™ 




Balmer/Paschen (NGC 7027) Jump 


11000 


7000 


8000 


8000 



a Versus [O II] temperature diagnostic. 
b Versus [N II] temperature diagnostic. 
c Versus [O III] temperature diagnostic. 
d Versus [CI IV] temperature diagnostic. 
c Versus [Ne III] temperature diagnostic. 



Table 4. Ionic Abundances a 



X+7H+ NGC 2440 IC 2501 IC 4191 NGC 7027 



Collisionally Excited Lines 



N°/H+ 


7 or+0.21 

' - zo -0.16 


6-9-6.6 


6-4-6.4 


6-1-0.3 


N+/H+ 


7 85+ ' 05 

' ■°°-0.06 


7 02+ - 16 


fi 74+ - 12 

°- '^-0.09 


7.19 


0°/H+ 


7 7fi +0.06 
'■'°-0.05 


7.68r6;i4 


7.34io.ii 


7.28_ 04 


0+/H+ 


7.24 


7.30 


7 c+0.3 
' -°-0.2 


7.68 


0+ 2 /H+ 


8-26lS:8S 


shifts 


8.79±0.05 


8.44 


Ne+ 2 /H+ 


7.4iiS:82 


7.58±0.03 


7 70+0.03 

' • ' y -0.04 




S+/H+ 


5.5±0.3 


5.93 


5.67 


5.80 


S+ 2 /H+ 


fi 23+ 007 


6.58lS;S 


6.08lg:S 


6.32±0.02 


C1+/H+ 




3.43 


3.18 


3.94 


Cl+ 2 /H+ 


4.70±0.06 


4-9218-.S 


4 56+ 007 


4.76±0.02 


Cl+ 3 /H+ 


/I fiO+0 06 
4 - 69 -0.07 


4.66t° if 


r 9 o+0.06 


4.72±0.02 


Ar+ 2 /H+ 


6.13±0.04 


6.27±0.07 


5 70 +0 ' 05 

°- ' u -0.06 


6.09±0.02 


Ar+ 3 /H+ 


5.70±0.03 


4.75±0.03 


r 74 +0.03 


5.77±0.02 


Ar+ 4 /H+ 
K+ 4 /H+ 


5.2418-.S 




4 ^7+0.09 

*-°'-0.12 


5.60 


4.59^;°? 




o 77+O.14 
° ' '-0.17 





Recombination 


Lines 








Hc+ /H+ 


10.86 


11.04 


11.00 


10.81 


Hc+ 2 /H+ 


10.72 


10.61 


10.08 


10.60 


C+ 2 /H+ 


8.65 


8.91 


8.46 


8.84 


N+ 2 /H+ 


8.43 


8.04 


8.28 


7.89 


0+ 2 /H+ 


8.51 


8.69 


9.08 


8.46 


Ne+ 2 /H+ 


8.08 


8.26 


8.76 





a In units of 12+log(X/H). 



Table 5. Atomic Data for Z>30 Ions 



Ion 



Levels 



Transition Probabilities 



Collision Strengths 



Br III ALL 2.05 

Biemont k Hansen (1986a) (2, 3) 

Br IV ALL 2.05 

Kr III ALL 2.05 

Kr IV ALL 2.05 

Kr V ALL 2.05 

Rb IV Persson k Wahlstrom (1985) 

Rb V Persson & Petterson (1984) 

Sr II Moore ('1952) 

Sr IV Hansen k Persson (1976) 

Sr V Ha nsen k Persson (1974) 

Sr VI Persson k Petterson (1984 ) 

Y V Reader k Epstein (1972) 

Y VI Persson k Reader (1986 ) 
Zr VII Reader k Acauista (1976 ) 

Te III Joshi. Tauheed. k Davison (1992 ) 

I III Tauheed k Joshi (1993b ) 

I V Kaufman. Sugar, k Joshi (1988 ) 

Xe III Persson et al. (1988) 

Xe IV Tauheed. Joshi. k Pinnington (1993 ) 

Xe V Gallard et al. (1999 ) 

Xe VI Chnrilov k Joshi (2000) 

Cs V Tauheed k Joshi (1993a ) 

Cs VI Tauheed. Joshi. k Kaufman (1991 ) 

Ba II Karlsson k Litzen (1999) 



Biemont k Hansen (1986a) 

Biemont k Hansen (1986a) 
Biemont k Hansen (1986b) 
Biemont k Hansen (1986a) 
Biemont k Hansen (1986a) 
Biemont k Hansen (1986b) 
Biemont k Hansen (1986a) 
Brage et al. (1998) (5-1, 5-2) 
Bie mont. Cowan, k Hansen (1988) 
Biemont k Hansen (1986b) 
Biemont k Hansen (1986a) 
Bie mont. Cowan, k Hansen (1988) 
Biemont k Hansen (1986b) 
Biemont k Hansen (1986b) 
Bie mont et al. (1995) 
Bie mont et al. (1995) 
Bie mont et al. (,1995) 
Bie mont et al. (1995) 
Bie mont et al. (1995) 
Bie mont et al. (.1995) 
Bie mont et al. (1995) 
Bie mont et al. (1995) 
Bie mont et al. (1995) 
Klose. Fuhr. k Weise (2002) 



Schoning (1997) 
Schoning (1997) 
Schoning (1997) 



Schoning k Butler (1998 ) 
Sch5ning k Butler (1998) 

Sch5ning k Butler (1998 ) 



Schoning k Butler (1998) 



Table 5 — Continued 



Ion 



Levels 



Transition Probabilities 



Collision Strengths 



;i993) 



Ba IV Sansonetti et al. 

Ba V Reader (1983) 

Ba VII Tauheed k Joshi (1992 ) 

Ba VIII Churilov. Joshi. fc Gavasov (2002) 

Pbll Moore ( 19581 



Biemont et al. (1995 ) 
Biemont et al. (1995 ) 
Biemont et al. (1995) 
Bie mont et al. (1995) 
Safronova. Safronova. & Johnson (2005) 



Schonine & Butler (1998) 
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Table 6. Emission Line Identifications and Intensities 



AT r (A) 


FWHMfkms" 1 ) 




F 

corr 


S/N 


Line Identification Notes a 








IC 2501 






3315.23 


23.73 


0.0088 


0.0200 


7.2 


[Mn III] 3315.01 :? 


3323.75 


19.16 


0.0116 


0.0260 


13.9 


Ne II 3323.73 


3327.17 


27.98 


0.0087 


0.0190 


10.4 


Ne II 3327.15 


3328.69 


19.13 


0.0056 


0.0130 


8.4 


N II 3328.72 


3334.83 


18.57 


0.0529 


0.1200 


50.0 


Ne II 3334.84 


3340.81 


13.75 


0.0188 


0.0420 


13.3 


N II 3340.87 


3342.64 


45.56 


0.0527 


0.1200 


37.0 


[Ne III] 3342.54 


3353.22 


33.40 


0.0071 


0.0160 


11.5 


[CI III] 3353.21 


7507.48 


22.85 


0.0048 


0.0027 


13.6 


N I 7507.61 : 


7509.97 


57.84 


0.0033 


0.0019 


7.1 




7512.69 


52.14 


0.0044 


0.0025 


7.1 


He I 7513.33 


7514.36 


36.15 


0.0038 


0.0022 


10.7 


C III 7514.30 


7515.77 


22.81 


0.0173 


0.0099 


32.2 


III 7715.99 


7530.42 


29.59 


0.5710 


0.3300 


627.6 


[CI IV] 7530.80 


7534.70 


30.21 


0.0022 


0.0013 


7.1 


[Xe IV] 7535.44 



Note. — The complete version of this table is in the electronic edition of the Journal. 
The printed edition contains only a sample. 

a : = uncertain identification; ? = uncertain feature; bl = line blend; ns = blend with 
night- sky line. 



Table 7. Possible Krypton Line Identifications 







NGC 2440 




IC 2501 


IC 4191 




NGC 7027 




IC 418 






Ao, 


5V a Mult°IDI a 


Ao, 


<SV Mult IDI 


A , SV Mult IDI 


A D , 


SV Mult IDI 


Ao, 


SV Mult IDI 




Transition 


I/Iff £!* 


(km s" 1 ) 




(km s^ 1 ) 


I/I H0 (kms- 1 ) 




(km s" 1 ) 


I/Ifffl 


(km s" 1 ) 


Ret 



6826.74 


1.8 


0/0 


5A 


6826.82 


5.3 


0/0 


5Abl 


6826.39? 


-13.6 


0/0 


6B: 


6826.91 


9.2 




. . -bl 


6826.87 


7.5 


0/0 


5Abl 




l.O(-B) 


-60.6 






4.6(-5) 


-57.1 






1.5(-5) 


-76.0 






5.7(-4) 


-53.2 






3.3(-4) 


-54.9 






B95.Z05 




-50.1 


0/0 




2.5(-5) 


-46.6 


0/0 






-65.5 


0/0 




4-6(-4) 


-42.6 


0/0 


■ • -bl 


S.2(-4) 


-AAA 






B95 




10.5 


4/0 


5A 




14.1 


4/0 


5A 




-4.8 


4/1 


3A 




18.0 


4/0 


7A 




16.3 


4/0 


7C 






-2.2 


4/0 


6C 




14.1 


4/0 


6C 




-17.6 


4/0 


9D 




5.3 


4/0 


7A 




3.5 


4/0 


6B 














-4.8 




■ • -bl 












12.7 




■ ■ -bl 








• • -bl 






-9.2 




... * 




























-1.8 












1 /l 


3A* 


5346 07 


2 8 


1 /1 


3A* 


004(1. UO 


2 2 


1 /l 


3 A* 


5345 99 


-1 7 


1 /1 


1A* 




-4 5 


i /i 


3A* 




3.5(-4) 


18.5 


2/0 


9 


6.1(-5) 


20.2 


2/0 


9 


2.4(-4J 


19.6 


2/0 


9 


1.9(-3) 


15.7 


2/0 


9 


3.5(-5) 


12.9 


2/0 


9 


H95,P04 




9.0 


5/0 


7 




10.7 


5/0 


6D 




10.1 


5/0 


6 




6.2 


5/0 


6 




3.4 


5/0 


5 


Z05 




-90.3 


8/0 






-88.7 


8/1 






-89.2 


8/0 






-93.2 








-96.0 






ZL02 




3.4 


0/0 


4B 




5.1 


0/0 


4B 




4.5 


0/0 


4B 




0.6 


0/0 


2B 




-2.2 


0/0 


5 






5.1 


/O 


4B 




6.7 


/o 


5C 




6.2 


/O 


5D 




2.2 


/O 


3C 




-0.6 


/O 


3A 






-11.8 


3/0 


7 




-10.1 


3/1 


6D 




-10.7 


3/2 


4B 




-14.6 


3/0 


8 




-17.4 


3/0 


8 




5867.68 


-3.1 


1/1 


3A* 


5867.82 


4.1 


1/1 


4A* 


5867.83 


4.6 


1/1 


4A* 


5867.71 


-1.5 


1/1 


1A* 


5867.73 


-0.5 


1/1 


2A* 




4.6(-4) 


2.0 






8.5(-5) 


2.0 






1.9(-4) 


2.6 






2.6(-3) 


3.6 


5/0 


5C 


3.5(-5) 


4.6 


4/0 


6 


B00.S03 




16.4 


6/0 


9 




20.5 


6/0 






21.0 


6/0 


>9 




14.8 


6/0 


8 




-13.3 


6/0 


8 


H01 




-15.8 








-8.7 








-8.2 








-14.3 


0/0 


7 




-13.3 


0/0 


7 


E04 




-68.5 








-61.4 


0/0 






-60.8 


0/0 






-67.0 








-66.0 






Z05 




36.3 


5/0 


7B 




43.5 


5/1 


0B 




44.0 


5/2 


4A 




37.8 


5/0 






38.9 


5/0 








0.0 








7.2 








7.7 








1.5 


0/0 


4B 




2.6 


0/0 


3B 




6107.68? 


-5.9 














6107.50 


-14.7 






6107.83 


1.5 


3/1 


5A* 












5.6(-6) 


19.6 


3/0 


6A: 


[4.3(-6)] 








1.9(-5) 


10.8 


3/0 


4A: 


7.2(-5) 


27.0 


3/0 




[1.4(-5)] 






































o /n 
1/0 


5A 














-13.8 
















-22.6 








-6.4 


9/0 


5A 












6798.01? 


-17.2 


2/0 


8 bl 


6798.15 


-11.0 


2/0 


7D 










6798.36 


-1.8 


3/1 


4Abl 












2.8(-5) 


-4.0 


7/2 


lAbl 


1.4(-5) 


2.2 


7/5 


1A* 


[8.1(-6)] 








4.1(-5) 


11.5 


7/1 


7bl 


[2.1(-5)] 








B95 


l.S(-5j 


-11.9 


0/0 


6 




-5.7 














1.6(-5) 


3.5 


0/0 


5B 












6256.31 


12.0 


0/0 


5Dbl 


6256.57 


24.5 


0/0 


7bl 


6256.31 


12.0 


0/0 


5D 


6256.38 


15.3 


1/0 


7bl 












4.4(-5) 


-41.7 


2/0 


8 


5.2(-5) 


-29.2 


2/0 


8 


1.5(-5) 


-41.7 


2/0 


8 


1.7(-4) 


-38.4 




■ ■ -bl 


[8.4(-5)] 








S03 


1.5(-5j 


-10.1 


2/0 


5D 


8.0(-6) 


2.4 


2/1 


2Abl 




-10.1 


2/0 


5D* 


1.2U) 


-6.7 


2/0 


6 bl 










Z05 




-1.9 


1/0 


4A 




10.5 


1/0 


6 




-1.9 


1/0 


4A 




1.4 


1/0 


5D 














-2.4 


5/0 


4A 




10.1 


5/0 


5D 




-2.4 


5/0 


4A 




1.0 


5/0 


4A 














-7.7 


4/0 


4A 




4.8 


4/0 


5D 




-7.7 


4/0 


5D 




-4.3 


4/0 


5D 














3.4 


4/0 


5D 




15.8 


4/0 


7 




3.4 


4/0 


4A 




6.7 


4/0 


7 














-6.7 


0/0 


5D 




5.8 








-6.7 


0/0 


5D 




-3.3 


0/0 


4A 














-7.2 


1/0 


5D 




5.3 


1/0 


4B 




-7.2 


1/0 


5D 




-3.8 


1/0 


4A 






















5.3 




■ • -bl 












-34.0 




■ ■ -bl 














-12.4 




■ ■ bl 




































OUT 








OUT 








OUT 








8244.33? 


34.2 






8243.70 


11.3 
































1.4(-4) 


70.9 






4.2(-4) 


47.7 






B95 




























23.3 


0/0 


5A 




0.4 


0/0 


2A* 


S03 



Notes 



[Kr III] 4p 4 3 P 2 -4p 4 1 D 2 A6826.70 

C I 3p 1 P i -4d i D 2 (V21) A6828.12 

Ho I 3s 3 Si-16p 3 P° A6827.88 

[Fc IV] 3d 5 4 P 1/2 -3d 5 2 D5 3/2 A6826.50 

Pc II] c4F 7/2 -y2G° /2 A6826.79 

SKY: OH 7-2 Ri(3.5) A6827.46 

SKY: OH 7-2 R,(4.5) A6828.47 



[Kr IV] 4p 3 



.3/2 



4p 3 



D ° /2 A5346.02 



S II 4s' ^D 3/2 -4p' ^F° /2 (V38) A5345.71 
C III 2s4d 3 D 2 -2p3d 3 Pf (V13.01) A5345..' 
[Fc II] a4F 3/2 -b4P 3/2 A5347.65 
Mg I 4s 3 Si-9p 3 P° A5345.98 
Pc II f4D 7/2 -4f 2[3]° /2 A5345.95 
Si II 6p 2 P 3/2 -4p' 2 P 3/2 A5346.25 



3/2 



3/2 
-4p 3 



[Kr IV] 4p 

Al II 4d °D-6f 3 F° (V41) A5867.64,.78,.89 
Si II 4s' 4 P° /2 -4p' 4 Pi/ 2 (V48) A5867.42 
Ni II 4f 2[4]° /2 -7d 4 H 9/2 A5867.99 
Ho II 5-29 5869.02 
[Cr III] 3d 4 5 D 3 -3d 4 3 H 5 A5866.97 
Ni II 4f 2[5]° 1/2 -7d 2 G 9/2 A5867.68 

[Kr IV] 4p 3 2 Dg /2 -4p 3 
[Fc II] a6D 9/2 -a2G 7/2 A6107.28 
Ca I 4s7d 3 D 2 -3d20p 3 P^ A6107.84 



[Kr IV] 4p 3 2 D 



C II 3s' *P° /2 -3p' *D 1/2 (V14) A6798.10 
Ca I 3d5s ^a-SdlSp A6798.28 



[Kr V] 4p 2 



-4p^ 



'D 2 A6256.06 



C II 4p ^P° /2 --> d ZD 3/2 (V10.03) A6257.18 

C II 3d' 2 D° /2 - 4p' 2 P 1/2 (V38.03) A6256.52 

S II 5p 4 S° /2 -7s 4 P 5/2 A6256.35 

Pc I a3H 4 -z3GJ (V169) A6256.36 

O I 3p' 3 F 4 -4d' 3 G I (V50) A6256.47 

C II 4p' 4 D 5/2 -5d' 4 D° /2 A6256.24 

Ca I 4s9s 1 S -3d25p 1 PJ A6256.45 

[Co II] a3D 2 -alF 3 A6256.46 

SKY: OH 9-3 Q 2 (0.5) A6256.94 

SKY: OH 9-3 Qi(1.5) A6257.96 

[Kr V] 4p 2 3 P 2 -4p 2 1 U 2 A8243.39 

___ T ___ 

H I 3-43 A8243.69 



Table 7 — Continued 





NGC 2440 


IC 2501 


IC 4191 


NGC 7027 


IC 418 


Transition 


X a , SV" Mult°[DI d 
I/I H fl"(k'» s" 1 ) 


Xa, SV Mult IDI 
I/lBf, (km s- 1 ) 


X a , SV Mult IDI 
I/I/tp (km s- 1 ) 


X„, SV Mult IDI 

(km s- 1 ) 


Xo, SV Mult IDI 
I/I H( 3 (km s _1 ) Ref Notes 


O III 5g G 2 [9/2]°-6h H 2 [11/2] A8244.10 








8.4 0/0 7B 


-14.6 B95 



References — B95 = PN NGC 7027 l lBaluteau et alJll995h. B00=Orion Nebu la {Baldwin et IDl2000F) . E04=Orion Nebula jEsteban et al~ll2004T) . H95 = PN NGC 6886 iHvung et alJll99St) . H01 = PN IC 5217 
iHyung et "alii 20011) , P04=PN NGC 5315 jPeimbert et al.||2004D , S03=PN IC 418 iSharpee et al.||2003l )'~ZL02 = PN Mz 3 iZhang &c Liu||2002fl , and Z05=PN NGC 7027 jZhang et alj|2005l) 

a Wavclcngth: 1) A D arc geocentric corrected wavelength in A, 2) "OUT" means not in observed range, 3) "?" denotes an uncertain feature. Intensity: 1) numbers in parentheses are exponents, 2) italics denote 
corrected intensities attributable to the s-proccss transition, 3) bracketed values as upper limits for unobserved features. 

^Observed (A )-transition wavelength (km s — 1 ). 

C EMILI multiplct check statistics: number expected/number observed 

EMILI IDI value/rank followed by asterisk (certain ID), colon (uncertain ID), or "bl" (blend). Definition of IDI given in Section R)l 
c Unidentified line in IC 418 |Baldwin et all feoOOl) . 
f Unidentified line in Orion Nebula fSharpee et all i2003|) . 
S NGC 2440, IC 4191, NGC 7027: Identified as a separate line. 
h NGC 7027: Affected strongly by telluric absorption. 



Table 8. Possible Xenon Line Identifications 







NGC 2440 






IC 2501 








IC 4191 








NGC 7027 






IC 418 










Ao, 


SV D 


Mult 


IDI d 


Ao, 


SV 


Mult 


IDI 


Ao, 


SV 


Mult 


IDI 


Ao, 


SV 


Mult 


IDI 


Ao, 


SV 


Mult 


IDI 




Transition 


T /J a 


(km s _1 ) 






T /J 


(km s ) 






T 11 

l/ljr/3 


(km s ) 






T /J 


(km s" 1 ) 






T It 

l/ljf/3 


(km s ) 






rtet IN otcs 


[Xe III] 5p 4 3 P 2 -5p 4 *D 2 A5846.77 


5846.73 


-2.1 






5846.64? 


-6.7 






5846.77 


0.0 






5846.65 


-6.2 




■ • -bl 


5846.70 


-3.6 


0/0 


5 * 


8 


Ho II 5-31 A5846.66 


3.2(-4) 


3.6 


0/0 


BB* 


7.7(-6) 


-1.0 


()/() 


3A 


5.2(-5) 


5.6 


0/0 


4A* 


4.2(-4) 


-0.5 


0/0 


lAbl 


1.3(-4) 


2.1 


0/0 


3A 


Z05 


[Fc II] a2G 7 / 2 -c2G 9/2 A5847.32 




-30.3 


3/0 


7C 




-34.9 


3/0 


7C 




-28.2 


3/0 


7G 


l-4(-4) 


-34.4 


3/0 






-31.8 


3/0 






Fc II 4d e4F 5/2 -4f 2[3]g, 2 A5846.78 




-2.6 


0/0 


3A 




-7.2 


()/() 


4B 




-0.5 


0/0 


4A 




-6.7 


0/0 


4B 




-4.1 


0/0 


4C 




rC 1 XOr — DQ ZyO/Z\\ A0040.DU 




6.7 
























2.6 










0/0 


3 A 




[Xe IV] 5p 3 4 S° /2 -5p 3 2 D° /2 A5709.20 


5708.95? 


-13.1 


()/() 


5 * 


5708.97? 


-12.1 


1/0 


5 : 


5708.98? 


-11.6 


1/0 


7* 


5708.89 


-16.2 


1/1 


7 * 


5708.91? 


-15.2 


1/1 


7* 


e 


N II 3s 3 Pg-3p 3 D 2 (V3) A5710.77 


1.2(-5) 


-95.6 






6.5(-6) 


-94.6 






1.2(-5) 


-94.1 






l.l(-4) 


-98.8 






1.9(-5) 


-97.7 






f 

many 


[Fc I] a5D 3 -a5Pi A5708.97 




-1.1 


7/0 


4C 




0.0 


7/1 


3A: 




0.5 


7/0 


4B 




-4.2 


7/0 


5C 




-3.2 


7/0 


5B 




Si I 4p 3 D 2 -18s (3/2,1/2)5" A5708.91 








3A 




3.2 


()/() 


3A 




3.7 


0/0 


3A 








3A 




0.0 


0/0 


4A 




Fo II] y4P° /2 -c6D 1/2 A5709.04 




4 7 


7 / 1 


3A 




-3.7 


7/0 


5 




-3.2 


7/0 


SD 




7 9 


7/1 
' / 1 


(iD 




-6.8 


7/0 


8 




[Xe IV] 5p 3 4 S° /9 -5p 3 2 D° /2 A7535.4 


7534.70 


-27.9 


0/0 


7 * 


7534.71 


-27.5 




■ ■ ■* 


7534.73 


-26.7 


1/0 


8 * 


7534.94 


-18.3 


1/1 


6 * 


7534.90 


-19.9 


1/1 


6 * 


g 


Fe II] b2Fgy 2 — z4Fg , 2 (V87) A7534.82 


1.3(-5) 


-4.8 


()/() 


2 A 


1.6(-5) 


-4.4 


1/0 


4 A 


3.6(-5) 


-3.6 


1/0 


3B 


1.7(-4) 


4.8 


5/0 


6 


1.6(-5) 


3.2 


5/0 


4 A 


H01 11 


N II 5f.G 2[7/2] 4 -10d 4 Fg A7535.10 




-15.9 


0/0 


6 




-15.5 


()/() 


6C 




-14.7 


0/0 


5 




-6.4 


0/0 


3A 




-8.0 


0/0 


5 


E04,P04 


No I 3p 2[l/2]i-3d 2[1/2]J (V8.01) A7535.77 




-42.6 


0/0 


7 




-42.2 








-41.4 








-33.0 








-34.6 








[Cr II] b4D 5/2 -c4D 5/2 A7534.80 




-4.0 


4/0 


5 




-3.6 


8/0 


6C 




-2.8 


8/0 


4C 




5.6 


8/0 


6 




4.0 


8/0 


4A 




No II 3d' 2 P 3/2 -6f 2[3]° /2 A7534.75 




-2.0 


0/0 


3B 




-1.6 


()/() 


BB 




-0.8 


0/0 


2A 




7.6 


0/0 


SD 




6.0 


0/0 


4A 




[Xe V] 5p 2 3 P -5p 2 3 P 2 A7076.8 


7077.02 


9.3 














7076.94? 


5.9 


0/0 


4A* 


7077.00 


8.5 


0/0 


6 : 












C I 3p 3 D 2 -4d J Dg (V26.01) A7076.48 


1.5(-5) 


22.9 






[5.5(-5)] 








3.9(-6) 


19.5 


6/0 


8 


2.0(-5) 


22.0 


6/0 


8 


[2.8(-5)] 








B95 


[Fc III] 3d 6 3 P4!-3d 6 1 S4 A7078.10 




-45.8 


0/0 


5A 












-49.2 


0/0 






-46.6 


0/0 












ZL02 


[Ni II] 4s 4 F 3/2 -4s 4 P 3/2 A7078.04 




-20.3 
















-46.6 








-44.1 














ZL02 


Ni II 5d 4 F 7/2 -4f 2[3]° , 2 A7076.90 




28.0 
















1.7 


0/0 


4A 




4.2 


0/0 


5C 












Ca I 4sl5s 3 Si-3d22p A7077.02 




22.9 
















-3.4 


0/0 


5C 




-0.8 


0/0 


4A 












Ca I 4sl4s 1 S -3dl7f L PJ A7077.08 




20.8 
















-5.5 


0/0 


5C 




-3.0 


0/0 


4A 












SKY: 3-2 2 6 1 E+-X 3 Eg Pi 2 (ll) A7078.58 




0.4 




. . .* 




































[Xe VI] 5p 2 P^,-5p 2 P 3/ / 2 A6408.9 


6408.69? 


-9.8 














6408.96? 


2.8 






6408.61 


-13.6 




... * 












Ho II 5-15 A6406.38 


2.8(-5) 


108.1 






[2.4(-6)] 








1.6(-6) 


120.8 






2.0(-4) 


104.4 






[9.3(-6)] 










[Fc III] 3d 6 1 S4 -3d 6 3 P2 2 A6408.50 




8.9 


1/0 


4A 












21.5 


1/0 


6A 




5.1 


1/0 


4A 












C IV 9-17 A6408.70 




-0.5 
















12.2 








-4.2 


0/0 


5B 













References. — E04=Orion Nebula iEsteban et al.ll2004H , H01 = PN IC 5217 iHvung et "aTll2001l) . P04=PN NGC 5315 jPeimbert et al.l 120041 ) , ZL02=PN Mz 3 jZhang fc Liul 120021 ). and Z05=PN NGC 7027 
(Zhang ot al. 2005) 

a ~ d Samc as for Tablc[7] 

"Unidentified lino in IC 418 jSharpoe et al.ll2003T) . 
1 Identified as separate line in all spectra. 

s On edge of NGC 2440, IC 2501, IC 4191 spectra, where poorest wavelength calibration is expected. 
h Listed bv lHvung et al.l I2001D as "unlikely ID". 
'IC 4191 and NGC 7027: Identified as separate lino. 



Table 9. Possible Bromine Line Identifications 







NGC 2440 




IC 2501 








IC 4191 








NGC 7027 






IC 418 






Transition 


A , 


SV" Mult c IDI d 

(km s _1 ) 


Ao, 


SV 
(km s _1 ) 


Mult 


IDI 


Ac 
yi-HfS 


SV 
(km s _1 ) 


Mult 


IDI 


Ao, 
I/I/f/3 


SV Mult 

(km s^ 1 ) 


IDI 


Ao, 

I/Ijf/3 


SV 
(km s _1 ) 


Mult 


IDI 

Ref Notes 


[Br III] 4p 3 4 S° / ,-4p 3 2 D° /2 A6131.0 






6130.70 


-14.7 


1/1 


4A* 


6130.34 


-32.3 


1/1 


6* 


6130.32 


-33.3 


■ ■ -bl 


6130.47 


-25.9 




. . .* <=■/ 


C III 7h ± '°H°-16g i o G A6130.30 

]Ni VI] 3d 5 4 D 5/2 -3d 5 2 Fl 7/2 A6130.40 


[1.7(-S)] 




1.7(-5) 


19.6 
14.7 


0/0 
1/0 


5B 
7 


3.8(-5) 


2.0 
-2.9 


()/() 
1/0 


2A 
5D 


l.l(-4) 
7.7(-S) 


1.0 0/0 
-3.9 1/0 


2Abl 
4C 


2.8(-5) 


8.3 
3.4 


()/() 
1/0 


2A 
4C 



[Br III] 4p 3 4 S 



Fc II 4d 4 P 



., 2 



3/2 



-4f 



O II 4f F 2 [4] 



9/2 



g 2 [5]ll/2 



^3/2 



AG555.94 
A6555.84 



6555.98? 
4.2(-5) 



N II 4f F 2 [5/2] 3 -6d 3 Dg A6556.06 
O II 4f F 2 [4]g /2 -6g 2 [5] 9/2 A6555.99 
O II 4f H 2 [5]° /2 11/2 -6g' 2 [6] AA6556.05,.08 
O II 4f F 2 [4]° /2 -6g 2 [4] 9/2 A6556.10 

[Br IV] 4p 2 3 P!-4p 2 'D 2 A7368.1 
C II 3p' 

O II 4p 1 y 2 

C V 7p 3 P°-8d 



-19.2 


1/0 


7 : 


6555.92? 


-22.0 


1/1 


6 * 


6555.87? 


-24.2 


1/1 


6 * 


6555.93? 


-21.5 


1/0 


8 * 


1.8 


0/0 


3 


7.3(-5) 


-0.9 


0/0 


2A 


8.5(-5) 


-3.2 


0/0 


2A 


4.3(-4) 


-0.5 


0/0 


2A 


(i.4 


0/0 


3 




3.7 


0/0 


2A 




1.4 


(1/(1 


2A 




4.1 


0/0 


3B 


-3.7 


0/0 


2 A 




-6.4 


0/0 


2A 




-8.7 


0/0 


3C 




-5.9 


0/0 


3D 


-0.5 


0/0 


2 A 




-3.2 


0/0 


3D 




-5.5 


(1/(1 


3C 




-2.7 


0/0 


3D 


-3.2 


0/0 


2A 




-5.9 


0/0 


3D 




-8.2 


0/0 


3C 




-5.5 


0/0 


4 


-5.5 


0/0 


2 A 




-8.2 


0/0 


4 




-10.5 


()/() 


4 




-7.8 


()/() 


5 



D 5/2 -3d' J P^ /2 A7370.00 
" -5s 2 P 3 /2 



p 3/2 A7367.6 
J D 7367.60 



[Br IV] 4p 2 3 P 2 -4p 2 1 D 2 A9450.5 

Fo I v3Dg-6d 2[5/2] 3 A9450.95 



[2.6(-S)] 



OUT 



[4.4(-5)] 



OUT 



[3.5(-5)] 



OUT 



7367.62? 
4.4(-5) 



9450.85? 
7.5(-5) 



-19.5 
-96.9 
-2.4 
0.8 

11.1 

-3.2 



1/0 9 



1/0 
0/0 



3 A 
4B 



1/0 8 : 
0/0 3A 



[2.1(-4)] 



[4.9(-5)] 



[5.3(-5)] 



References. — B95 = PN NGC 7027 iBaluteau et al.lll995l l 
a - d Samc as for Tablc[7] 

c Computed using energy levels from Bicmont Hansen ( 198 6al) . see text. 
f Unidentified line in IC 418 lSharpee et al.l 12003D . 



Table 10. Possible Identifications From Other Z> 30 Ions 







NGC 2440 




IC 2501 






IC 4191 








NGC 7027 






IC 418 












SV D 


Mult' 


IDI a A„, 


SV 


Mult IDI 




SV 


Mult 


IDI 




SV 


Mult 


IDI 


Ao, 


SV 


Mult 


IDI 




Transition 


T /T ^ a 


' kmS - 




I/I/f/3 


S - 




I/I/f/3 








I/I/f /3 








I/I/f/3 


(kms } 






JTvCI 1NOTCS 


[Rb IV] 4p 4 3 P 2 — 4p 4 1 D2 A5759.55 


5759.98 


22.4 










5759.61 


3.1 




■ • -bl 


5759.59 


2.1 




• - -bl 


5759.78? 


12.0 






e,/ 


He II 5-47 A5759.74 


5.1(-4) 




0/0 


6B* [2.5(-6)] 






2.4(-5) 


-6.8 


0/0 


4Bbl 


3.1(-4) 


-7.8 


0/0 


5Abl 


2.0(-5) 


2.1 


()/() 


3A 


Z05 


[Fc II] a2H 9/2 -c2D 5/2 A5759.30 




3j.4 


0/0 


5A 






1.5(-5) 


16.1 


0/0 


3A 


z.0(-4) 


15.1 


0/0 


6B 




25.0 


()/() 


5B 




[Rb V] 4p 3 4 S° /2 -4p 3 2 D° /2 A5363.6 








5363.81 


11.7 




5363.96? 


20.1 






5363.62 


1.1 




. . .* 












[Ni IV] 3d' ^F 7 / 2 — 3d' ^Gg/ 2 A5363.35 


[3.4(-5)] 






6.7(-6) 


25.7 


4/0 7 


5.4(-6) 


34.1 


4/0 


8C 


9.0(-5) 


15.1 


4/0 


7B 


[2.3(-5)] 








many 


O II 4f.F 2 [4]° /2 -4d' 2 F 7/2 A5363.80 










0.6 


0/0 3A 




8.9 


0/0 


5A 




-10.1 


0/0 


6A 
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Table 11. Summary of Possible Z>30 Ion Line Intensities' 1 
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NGC 7027 



Transition 
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[Br TTTl Afi1 31 f) 




1.7(-5) 
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4.1( 
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<1.4(-5) 






[Kr IV] A6798.4 
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Baldwin et all (|2000h : NGC 3576= lGarria-Roias et al.l (|2004f ) 



''In units of I(H/3)=1, certain identifications in bold. 
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Table 12. s-process Ionic/Elemental Abundances' 









PNe 






H II regions 


X+7H+ /Transition 


NGC 2440 


IC 2501 


IC 4191 


NGC 7027 


IC 418 


Orion 


NGC 3576 


Br+ 2 /H+ A6131.0 
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9 91 +0.16 
z - zi -0.19 


9 -, o+0. 14 
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9 r 9 +0.09 
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9 1 O+0.19 
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9 1 O+0.19 
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A5867.74 




2.fl9±g : g 
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2 4 n+o.i6 
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q O+0.3 
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<2.60 
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Ba +3 /H+ A5696.6 
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1 7+0.3 
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6.13±0.04 
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c 7n +0.05 
°- '"-0.06 


6.09±0.02 


fi 01 + - 06 
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6.42 


6.34±0.04 b 
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5.70±0.03 


4.75±0.03 
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5.77±0.02 
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z -yo_o.o6 


4.37 


4.20±0.07 b 
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°- z *-0.08 




4 C7+0.09 

4 - ' -0.12 


5.60 








Br/Ar 
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Kr/Ar 


-3.33 ±0.11 


-3.3 ±0.3 


_ 2 Q4+0-09 
z - y ^-0.10 


-2.36 ±0.04 


-3.19 


< -3.50 


Xe/Ar 


-4 64+ 012 


_ 4 01+0.15 

4 - oi -0.19 


O 7O + 0.13 


q qr+0.08 
°- Jd -0.09 


_o 04+O.W 


< -3.29 ••• 


[Br/Ar] 








,n 70+0. 16 
i-u./o_ 22 








[Kr/Ar] 


-0.07 ±0.11 


-0.1 ±0.3 


+0 32+ ' 09 

+UiW -0 .10 


+0.90 ±0.04 +0.82+H 4 


+0.07 


< -0.24 


[Xe/Ar] 


-0 T7+ - 12 

U - J '-0.14 


-o.o4t8-.ig 


+0.49l^ 3 


+0 92+ - 08 

' u - az -0.09 


+o.93t8;if 


< +0.98 • • • 



a In units of 12+log(N/H). 

b Ar abundances using t 2 =0.00. 
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Fig. I. — Portion of the IC 2501 spectrum illustrating the automated continuum fitting 
procedure. The vertical lines designate emission features defined by the ED fitting algorithm. 
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Fig. 2. — The fraction of lines present in the spectrum of each PN which are also detected 
in another of the PNe. The numbers give the total number of detected lines in each flux 
interval. 
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Fig. 3. — The cumulative number of observed lines exceeding a given flux level for recently 
published nebular spectra. We consider only those lines within the wavelength range 3510- 
7470 A, which is covered by all the spectra. 
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Fig. 4. — Continuum subtracted spectra in the immediate vicinity of several prominent Kr, 
Xe, and Br lines, shifted to each nebula's rest frame. Predicted wavelengths are indicated 
by dashed lines. Labels "NS" indicate telluric night glow lines mentioned in the text. The 
regions around [Kr III] A6826.70, [Kr V] A6256.06, and [Xe III] A5846.77 are expanded in 
subsequent figures and discussed in the text. 
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Fig. 4. — Continued. 
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Fig. 4. — Continued. 
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Fig. 4. — Continued. 
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Fig. 5. — LEFT: Spectra in the vicinity of [Kr III] A6826.70 (black line) shifted to nebular rest 
wavelengths, with similarly shifted and superimposed coadded modeled spectra of telluric 
emission and observed interloping nebular features (dot-dashed gray line). Pluses indicate 
telluric OH Meinel 7-2 Rx(3.5) A6827.46, Ri(4.5) A6828.47, and Ri(2.5) A6829.49 lines (left 
to right), nebular He I A6827.88 is marked with an arrow, and a O VI A1032 Raman feature 
at A6829.16 (NGC 7027) marked with a lower arrow. RIGHT: Same for [Kr V] A6256.06, 
with OH Meinel 9-3 Q 2 (0.5) A6256.94 and Qi(1.5) A6257.96 marked as pluses (left to right), 
and nebular C II A6256.52, A6257.18, A6259.65 lines marked with arrows (left to right). 
Panels to right show residuals after subtraction in vicinity of predicted wavelengths of target 
line (vertical dashed line in all panels). 
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Fig. 6. — Spectra in vicinity of [Xe III] A5846.77 (black line) shifted to nebular rest wave- 
lengths, with superimposed profile of the He II 5-32 A5837.06 shifted to the rest wavelength 
of the He II 5-31 line and scaled as described in the text (gray dot-dashed line). Panels to 
right show residual spectra after subtraction near the rest wavelength of the [Xe III] A5846.77 
line (vertical dashed line in all panels). 



